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In part I general aspects of the renormalization of a spontaneously broken gauge theory have 
been introduced. Here, in part II, two-loop renormalization is introduced and discussed within 
the context of the minimal Standard Model. Therefore, this paper deals with the transition 
between bare parameters and fields to renormalized ones. The full list of one- and two-loop 
counterterms is shown and it is proven that, by a suitable extension of the formalism already 
introduced at the one-loop level, two-point functions suffice in renormalizing the model. The 
problem of overlapping ultraviolet divergencies is analyzed and it is shown that all counterterms 
are local and of polynomial nature. The original program of 't Hooft and Veltman is at work. 
Finite parts arc written in a way that allows for a fast and reliable numerical integration with 
all collincar logarithms extracted analytically. Finite renormalization, the transition between 
renormalized parameters and physical (pseudo-)observables, will be discussed in part III where 
numerical results, e.g. for the complex poles of the unstable gauge bosons, will be shown. An 
attempt will be made to define the running of the electromagnetic coupling constant at the 
two-loop level. 
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1 Introduction 



After the end of the Lep era it became evident that including estimates of higher-order radiative cor- 
rections into one-loop calculations for physical (pseudo-)observables could not satisfy anymore the need for 
precision required by the forthcoming generation of experiments. Since LHC will be an arena for discovery 
physics, high precision will not be mandatory during its first phase. However, according to some predesti- 
nate design, hadron machines are alternating with electron-positron ones and, hopefully, ILC will come into 
operation: here the highest available theoretical precision will play a fundamental role. 

As a matter of fact, it is not clear what kind of scenario will arise after the first few months of running at 
LHC. Any evidence of new physics will lead to a striking search for new theoretical models and their Born 
predictions; the quake could be so strong to remove any interest in quantum effects of the Standard Model. 
On the contrary, we could be back to the familiar landscape: effects of new physics hidden inside loops. 

Since we have no firm opinion, we decided to follow the old rule si vis pacem para bellum, building the 
environment that allows for a complete two-loop analysis of a spontaneously broken gauge field theory. This 
construction requires several elements, and it is difficult to characterize our strategy with a single acronym; 
although our work implies a lot of analytical aspects, the final step (computing arbitrary two- loop diagrams) 
can only be done with an algebraic-numerical approach. 

If one thinks for a while, everything is in the old papers of 't Hooft and Veltman [1 , but translating few 
formal properties into a working scheme is far from trivial. Most of the times it is not a question of How do 
I do it?, rather it is a question of bookkeeping: Can I do it without exhausting the memory of my computer?, 
or, Is there any practical way of presenting my results besides making my codes public?. 

We devoted a first paper [2] (hereafter I) to deal with general aspects of a spontaneously broken gauge 
theory. First of all, we showed how to treat tadpoles; although everybody knows how to do it, general results 
are rarely presented in a way that everyone can use them. In addition, we analyzed how to perform an 
order-by-order diagonalization of the neutral sector of a theory of fundamental interactions; once again, one 
needs a comprehensive collection of results which allows for practical applications. 

Alternative solutions to solving these problems exist, noticeably in the background- field formalism [3] 
(compare also with Ref. [4 ) ; our claim here is restricted to the construction of a set of procedures which do 
not rely on other sources and cover broadly all the aspects, from generation of diagrams and renormalization 
to evaluation (mostly numerical) of physical processes. In particular, our results show that the structure of 
the counterterms at the two-loop level, as well as of the whole set of renormalized Green functions and the 
Ward-Slavnov- Taylor identities [5] , has (in the conventional approach) a degree of simplicity comparable to 
the one obtained at one loop in the background-field approach. 

Another perennial question is: what about renormalization, with or without counterterms? In a way, it 
is a fake question, since the two approaches are fully equivalent as far as S'-matrix elements are concerned. 
In this paper we focus on the transition from bare parameters to renormalized ones, and in a third paper 
(hereafter III) we will then discuss the ultimate step in any renormalization procedure: the transition from 
renormalized parameters to a set of physical (pseudo-)observables. Perhaps, one should try to make a 
clear vocabulary of renormalization in quantum field theory. A renormalization procedure is designed to 
bring you from a Lagrangian to theoretical predictions; it includes regularization (nowadays dimensional 
regularization pQ is easy to understand), a renormalization scheme and the choice of an input parameter 
set. The scheme, being a transitory step, is almost irrelevant; it can be on-mass-shell or MS or based on 
complex poles, but unless you do something illegal (resummations that are not allowed or similar things) it 
really does not matter. Admittedly, one can define MS quantities as convenient tools, but it is the last step 
that matters, at least as long as we have a convenient subtraction point (which we miss in QCD). 

Renormalized quantities should always be expressed in terms of a set of physical quantities. One may 
indulge to the introduction of an MS running electromagnetic coupling constant (importing a concept from 
QCD to QED, which sounds strange anyway) but, at the end of the day, only cross sections matter. 

In this paper, we have done one thing: all the two-loop Green's functions of the theory are made finite 
by introducing non-logarithmic counterterms and respecting unitarity. In addition, one can easily check that 
renormalized Ward-Slavnov- Taylor identities are satisfied. Actually, we have done more, since all ultraviolet- 
finite parts have been classified and an algorithm has been designed for their evaluation at any scale. What 



1 



is innovative in our approach, as well as in other modern approaches [5J, is the idea that everything can be 
generated (is generated) by a set of automatized procedures which deals satisfactorily with the somewhat 
greater complexity of a two-loop calculation. 

Furthermore, classification of ultraviolet divergencies is dynamically linked to a well-defined computa- 
tional scheme. In other words, in our approach, the ultraviolet-finite parts are back in their privileged 
position where they can play the role of creating the predictive power of the theory. There are, of course, 
preliminary steps (not always the easy ones), but it is only the full control on the multi-scale level that pays 
off. 

The outline of this paper is as follows. After introducing our notation and conventions in Section [21 
we outline the strategy of our calculation in Section [3l where we classify all the ultraviolet-divergent parts 
of the needed loop integrals. Next, in Section [5J we choose the gauge-fixing Lagrangian and we define our 
renormalization scheme. In Section [5] we review one- loop renormalization, and in Section [BJ we show explicit 
results for tadpole renormalization and neutral-sector diagonalization at two loops. Finally, we analyze 
the ultraviolet structure of two- loop self-energies in Section [7] and we show explicit results for the two- loop 
counterterms in Section [H Section [9] contains the summary, and several technical details connected with the 
relevant kinematical limits, which represent the backbone of our computational techniques, are discussed in 
the appendices. 



2 Notation and conventions 

Regularization. We employ dimensional regularization [I] , denoting the number of the space-time dimensions 
by n = 4 — e. In addition, we use a short-hand notation for regularization-dependent factors, 

M 2 M 2 
A uv — 7 + ln7r + In— 5-, A uv (x) = A uv - In , (1) 

fi z x 

where 7 = 0.5772156 • • • is the Euler constant, \i is the 't Hooft unit of mass, M stands for the bare 
(renormalized) VF-boson mass (we do not distinguish unless strictly needed) and a; is a positive-definite 
kinematical variable. In our conventions the logarithm has a cut along the negative real axis and it is 
understood that for all masses: M 2 — > M 2 — id, with 6 — > 0+. 

Masses. We introduce a compact notation for ratios of squared masses, 

_ M 2 H _ mf ti _ _ 

x H - A/2 , xi ti - M2 , x u>i - M2 , x d ,i - Af2 ■ V) 

Here M H is the Higgs-boson bare (renormalized) mass and nij <, m Uy i and m^i are the bare (renormalized) 
masses of the charged lepton and the up and down quarks of the ith fermion doublet, with i = 1, ... ,3. 
We consider the minimal representation for the Standard-Model scalar sector, defining Mq = M/cg for the 
Z-boson bare (renormalized) mass. Here eg (sg) is the bare or renormalized cosine (sine) of the weak- mixing 
angle 6. For notational clarity we frequently employ the notation 

{m}i2..jv =mi,m 2 ,...,m N . (3) 



3 Outline of the calculation 

Our calculation builds upon an automatic strategy for generating Feynman diagrams and evaluating the 
necessary one- and two-loop integrals. In addition, it does not rely on any black-box tool: diagrams are 
generated through a set of FORM [7] routines implemented in the QraphShot package [8] and loop integrals 
are computed through the FORTRAN/95 CoopBack [9] code. 

The present work uses a set of results derived in a series of previous papers. The general strategy for 
handling multi-loop multi-leg Feynman diagrams was designed in Ref. |10j and the whole body of results 
necessary for evaluating two- loop two-point integrals can be found in Ref. [llj . The calculation of two- loop 
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three-point scalar integrals is considerably more involved: infrared-convergent configurations are discussed 
in Ref. [12] and infrared- and collinear-divergent ones are analyzed in Ref. [13j . Finally, our method for 
dealing with two-loop tensor integrals can be found in Ref. [T3] and results for one-loop multi-leg integrals 
are shown in Ref. [T5] , 

It is worth noting that in any spontaneously-broken gauge theory the large number of masses and external 
kinematical variables leads to difficulties for the familiar analytical methods. Therefore, we use an alternative 
approach (mostly numerical) which, however, does not imply a blind rejecting of analytical techniques; for 
instance, in this work, we devote a special attention to the analytical extraction of leading and sub-leading 
collinear logarithms, since they play a dominant role in any calculation [16 . 

If we focus on renormalization, we can safely state that all the necessary ingredients are available. Here 
the crucial point is to connect a set of input experimental data (an input-parameter set, hereafter IPS) to 
the free parameters of the theory: 

- mass renormalization involves the calculation of self-energies jlOll lj ; 

- renormalization of coupling constants requires additional elements, which depend on the choice of the 
(pseudo-)observables (hereafter POs) in the IPS. 

The most-obvious selection of an IPS is based on the choice of those data which are known with the best 
experimental precision: the electromagnetic coupling constant, a, and the Fermi coupling constant, G F . 
^From the diagrammatic point of view, the whole set of renormalization equations (including finite parts) 
requires the evaluation of two-loop vertex and box integrals, since a is defined by means of the Thomson- 
scattering amplitude and G F is extracted from the muon lifetime. However, it is possible to show that: 

- neglecting O (rn^/M 2 ) terms, where m M is the muon mass, the evaluation of the muon lifetime can be 
carried on through tadpole diagrams (zero external momenta); 

- the relevant Ward-Slavnov- Taylor (hereafter WST) identities of the Standard Model imply a set of 
algebraic cancellations involving vertex diagrams and wave-function renormalization factors for the 
external legs. Here the role played by the parameter V introduced in I (all-order diagonalization of the 
neutral sector) plays an essential role. As a result, one can write an electric-charge renormalization 
equation which depends only on the photon vacuum-polarization function, as in QED. 

We defer a complete discussion of this computation to a forthcoming paper (hereafter III); for a short review, 
see Ref. [17]. For the present purposes, it is enough to say that the results collected in Refs. [101 1 lj allow to 
complete the whole renormalization program. Of course, properties of two-loop vertex diagrams have to be 
checked, and here we use the results of Refs. [12|13ll4) . 

An important remark concerns tensor loop integrals, namely those configurations with a non-trivial spin 
structure for the numerators. On the one hand, they can be reduced to scalar integrals if one uses the 
standard methods of Ref. [18] , supplemented by the sub-loop techniques of Ref. [19] . It is worth noting that 
a complete scalarization for two-loop vertex diagrams requires the introduction of generalized scalar integrals, 
to be evaluated in shifted space-time dimensions (the corresponding propagators are raised to non-canonical 
powers). On the other hand, the number of terms generated by the reduction procedure increases rapidly in 
any realistic computation. 

In order to carry out a systematic check of our computation, we have used both approaches: a reduction 
procedure followed by an extraction of the ultraviolet (hereafter UV) resides of the poles for the scalar 
configurations and a direct computation of the UV divergencies for the tensor cases. Reduction techniques 
can introduce negative powers of Gram determinants, which can be dangerous for the numerical stability 
of the results. The approach we take in the present paper avoids these problems since scalarization is only 
employed to prove that the relevant algebraic properties are satisfied. 

One of the main goals of this paper is to show that UV divergencies at two loops can be removed through 
suitable polynomial subtraction terms. Therefore, in this section we collect the results for the residues of 
the UV poles of all the relevant one- and two-loop integrals employed in our calculation. However, we stress 
that only a full control over the (less-trivial) UV-finite parts will finally allow for theoretical predictions. In 
Subsection 13. II we introduce general UV decompositions for one- and two-loop integrals. In Subsection 13.21 
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we present the relevant results for the one-loop ingredients. In Subsections 13.31 13.41 and 13.51 we show the 
necessary results for one-, two- and three-point two-loop integrals. 



3.1 Ultraviolet decompositions 

One-loop integrals. Any one-loop integral f 1 can be decomposed as0 



fH{i})= E /'(W; *) ( 4 ) 

fc=-i 

where {1} denotes a given set of arguments: powers of the inverse propagators, external kinematical variables, 
masses of internal particles, x is some kinematical variable (usually M for tadpoles or a squared external 
momentum) and the dependence on the dimensional regulator e and the regularization-dependent factors 
introduced in Eq.|T]) is entirely transferred to the universal UV factors, 

F^(x) = ^-^A uv (x) + ^A 2 uv (x)e, F£(x) = 1 - \ A uv (x) e, F}(x)=e. (5) 

It is worth noting that, because of overlapping divergencies (UV-divergent one-loop sub-diagrams), we include 
O (e) terms in all one- loop results. 

Two-loop integrals. A generic two-loop integral f 2 can be written as 

f({l})= E f 2 ({l};k)F 2 (x). (6) 

fe=-2 

Here the two-loop UV factors read as follows: 



F 2 2 (x) = ^-^^- + ^A 2 uv (x), Fl 1 (x) = - £ -A ov (x), F l7 '(,)=l. (7) 



Note that the product of two one-loop integrals can be written through the same UV decomposition of a 
two-loop integral. Finally, in the following we collect two overall factors for one- and two-loop integrals, 

Mi(e) = ^, M2 (e) = ^ = - M ?(e). (8) 

3.2 One-loop integrals 
Tadpoles 

A scalar one-loop tadpole, with the inverse propagator raised to an integer power a, is defined by 

d n q 
(q 2 + m 2 ) c 



A)(a,m) = m{e) / . (9) 



Integration-by-parts (IBP) identities [20] allow to reduce any configuration with a > 1 to the case a = 1 by 
means of a recursive application of the identity 

. . l/Ti \ Ar)(a — l.m) 

A (a,m) = + 1 — — ( 10 ) 



a - 1 V2 

We factorize everywhere the scale dependence in order to deal with dimensionless loop integrals. Here we 
introduce the function ao(a,m), 

A (a,m) = (m 2 ) 2 ~ a a (a,m), (11) 



lr The UV decompositions of Ecjs . ((4j— (TJl and the corresponding classification of UV poles are due to S. Uccirati. 
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and apply the UV decomposition denned in Eq.(|3]) to the case ao(m) — ao(l,m), 



«o 



(m)= J2 «o(m; k)Fl{M 2 ). 



(12) 



k=-l 



The coefficients of the UV factors in Eq. (|T2jl read as 



a (m; -1) 
a (m; 1) 



a (m; 0) 
-k( 2 )+ (' 



-1 + ln 



In 



M 2 ' 
to 2 \ 



In- 



M 2 / M 2 



(13) 



where C(z) is the Riemann zeta function. 
Two-point integrals 

Ultraviolet decompositions. For two-point integrals we introduce s = —p 2 , where p is the external momentum. 
We extract form factors for tensor configurations and introduce scaled quantities according to the following 
definitions: 



B lllJ {a,f3, s,{m} 12 ) = ^i(e) / d n q 



(q 2 + m\) a [{q+p) 2 + m^f 



b 2 i{a,f3,s,{m} 12 )p f _ l p v 



b 22 (a,/3, s,{to}i 2 )<5 m „, 



B M (a,^, s, {m}i 2 ) = /Lti(e) / d"g 



_ 2-Q-0 



S (a,/3, s,{m} 12 ) = m(e) 



(q 2 + m 2 ) a [{q+p) 2 + m 2 ]P 
d n q 

(q 2 + m\) a [{q + p) 2 + m 2 .]' 3 



&i(a,^,s,{TO}i 2 )p M) 



2-Q-/3 



b (a,/3, s,{m}i 



(14) 



Note that scaled one-loop form factors will often appear in the residues of the UV poles for two-loop integrals. 
Reduction of higher-rank form factors to the scalar case can be obtained by means of well-established 
techniques. UV decompositions read as 



bi{a,P, s,{m}i 2 ) = ^2 s ' { m li2; k ) F k( s )> 



(15) 



fc=-i 



where i = 0, 1, 21, 22. The coefficients of the UV factors for the canonical scalar case (a — (3 — 1) read as 



6 (1, 1, s, {m} 12 ; -1) = 2, b (l, 1, s, {m} 12 ; 0) = fe[j n (l, 1, a, {m} 12 ), 
6 (1, 1, «, {m} 12 ; 1) = 1 [C(2) - 6g(l, 1, s, {m} 12 ) 



(16) 



We will need to consider additional cases, where one propagator is raised to a non-canonical power, i.e. 
a = 2,(3 = 1. Here results are given by 



b (2, 1, «, {m}i 2 ; -1) = 0, 6 (2, 1, fl, {m} 12 ; 0) = ^ n (2, 1, s, {m} 12 ), 
6 (2, 1, s, {m} 12 ; 1) = - i 6g(2, 1, s, {m} 12 ). 



(17) 



Note that Eq. (fT7|) holds strictly only when mi ^ 0, otherwise an infrared (IR) divergency arises and we need 
a different decomposition. The finite parts and the O (e) terms introduced in Eas. !!^ and (TIT)) have the 
following integral representations: 



C(l,l, S ,{m} 12 ) = 



cfa; In 



C(2,l,s,{m} 12 ) 



dx 



{l-x)s 
X(x) 



6q(1, l,s,{m}i 2 ) = - / cb ln ; 



, 2 xfc) 



/.,',(2.l..s.W,,) = 2 / dx {l ^ S ln^, (18) 



5 



where the quadratic form x reads as 

x{x) — s x 2 — (s + m\ — m\ ) x + m\ — id. (19) 

Finite components. For the finite parts we compute explicitly the integrals introduced in Eq. (|18[) . For the 
case a = j3 = 1 we obtain 



6^(1,1,5,0,0) = 2-L 



6« n (l,l,s,0,m) 
6§ n (l, 1, s, m, m) 

6o n (l, l,S,TOl,TO 2 ) 



an ) 
2 



In h 2 - 1 In 1 



, m2 r, n,™? x , (3(m 2 /s) + 1 

ln_ + 2 -/3 — In A ' -, 

s s p(m z I s) — 1 



? 9 2 9 

, tot , m% tot — mS 
In — - + In — - H - 

s s s 



fln^i 

V s 



-In ^ 

s s 



+ 2- -A 1/2 (s,m?,TO2)i( s ,m?,TO2), 
where A(x, y, z) is the usual Kallen function and we introduced 



(20) 



a 2 { 



rn 



1 -4- 



s 

,2 „,2\ 



L(s, TOj, m 2 ) = In 



s + i5 

■S + ml + 1712 — -^ 1/,2 ( s 7 m l; m \) 



2 TOl TO2 

L aa = ln(-l-i*), <5^0+. 
For the regular (mi 0) case, where a = 2 and = 1, we obtain 



6[j n (2,l,s,m,0) = -ln(l- 



6g n (2,l, s,m,m) 



1 



/3(to 2 /s) + 1 
/3(to 2 /s) /3(to 2 /s) - 1 : 



In 



m? — m| — s . 2 2 . 1 / to? 

, s, TOi, to 2 ) = -r^ s 2- L(s, to 1; to 2 ) + - I In 

A 1 /^(s, m\, m\) 2 V s 



In- 



(21) 



(22) 



Infrared- divergent configurations. Two special cases requires a separate discussion: when toi = in Eq. (|14p 
and a = 2 or (3 = 1, an IR divergency shows up. Here we regularize IR divergencies by setting n = 4 + e, 
and then we use the relation between the UV regulator e and the IR one e, i.e. e = — e. The results read as 
follows: 



6o(2,l,s,0,0) 
6 (2, l,s,0,m) 



2 . . . e 

Apv^J - i an + - 

e 4 



A 2 ,( s ) + 2A^( s )L an + L 2 n -C(2 



r z . . , , to / to \ , / s 

{ A uv (s) — In IH 1 ln(l 

I e s \ s I \ to^ 



- C(2) + - A£ V ( S ) + In — A uv (s) + - In — 



1 , to 



1 + — In 1-— A uv (s) + ln— hi--; 



rn- 



- In 2 1 - — - Li 2 



)]}■ 



Here we used generalized Nielsen polylogarithms, 

(_1)«+P-X yl ^ ^n-i 



(n- l)!p! y a; 



ln p (l-za;), S n _ 1,1(2) = Li„ (2) . 



(23) 



(24) 



G 



Three-point integrals 

Ultraviolet decompositions. For three-point integrals we extract form factors for tensor configurations and 
then we introduce scaled quantities. Our notation for a generic tensor structure with I Lorentz indices is 

3 W • • • %i 



C» 1 ... n (a 1 ,a 2 ,a a ,p 1 ,p 2 ,P,{m} 123 )=» 1 (e) J d» q + ^ )2 + ^ [fa + p)2 + ^ , 

where P = pi + p 2 . Scaled quantities read as (we omit list of arguments) 



Cuii 



_ 2-a 1 -a 2 -a 3 



C2lPlp.Plv + C22P2^iP2v + C 2 3 {P1P2} fiv 



= s'-<*i-<**-<*3 {c nPlfi + c 12 p 2fi ), C = s z -^- a ^c 0l (26) 

where we have also introduced the symmetrized combination {pk}^ — p^ k v + p v k^. Three-point, scaled, 
one- loop form factors will often appear in the residues of the UV poles for two-loop vertex integrals. Reduc- 
tion of higher-rank form factors to the scalar case can be obtained by means of standard methods, and UV 
decompositions read as 

1 

c i (a 1 ,a 2 ,a 3 ,p 1 ,p 2 ,P, {m}i 23 ) = ^ Cj(ai, a 2 , a 3 ,p 1: p 2 , P, {m}i 23 ; k) F£(s), (27) 

fc=-i 

where i contains the information on the spin structure and s = —P 2 . Examples of UV decomposition are 
c (l,l,l,pi,p 2 ,P,{m}i 23 ; -1) =0, ... c 24 (l,l,l,pi,p 2 ,P,{m}i 23 ; -1) = -, etc. (28) 

3.3 Two-loop tadpoles 

A special class of two-loop integrals is represented by tadpoles, defined by 

T ijk «m} 123 ) = M2 (e) / — , d "^ f g2 — . (29) 

J {qf+mf) {q 2 2 +ml) ] [{ qi - q 2 ) 2 + mlf 

Any two-loop tadpole can be reduced to a single two-loop master tadpole (we choose the configuration with 
i = k = 1, j = 2) and to products of one- loop tadpoles through iterated application of IBP identities. A few 
examples are given in Appendix fAJ The UV decomposition for the master tadpole reads as 

o 

T m ({m} 123 ) = T 121 ({m} 123 ;k) F 2 (M 2 ), (30) 

fe=-2 

where the coefficients of the UV-divergent factors are given by 

2 

7\ 2 i({m}i 23 ; -2) = -2, T 121 ({m} 123 ; -1) = - ( 1 - 2 In ^| ) . (31) 

These relations hold if m 2 ^ 0, since a collinear logarithm is present when m 2 — > 0. When m 2 = we can 
change integration variables and apply again IBP identities in order to use Eq.(j3Tjl. i.e. 



T121 (mi, 0,7713) = T112 (mi, 7713,0) — i • Tm (mi, 77i 3 ,0), (32) 



IBP 



and the (mi,m 3 ,0) configuration, which is collinear-free, can be computed as shown in Eq. (|223[ ) of Ap- 
pendix [A] 
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3.4 Two-loop two-point integrals 

Notation. Wc follow the notation of Section 7 of Ref. [2] and we introduce scaled irreducible two-loop 
two-point integrals. Note that when the squared external momentum vanishes, p 2 = 0, the tadpole-reduction 
methods of Subsection 13.31 can be easily employed. Otherwise, we extract s = —p 2 and accordingly define 



SSq 4 , 



D 



b > 



CD 

°0 



1 



QE 

D 



1 



The alphameric classification of graphs can be found in Subsection 2.2 of Ref. 
we employ form factors according to Section 7 of Ref. [14] , 



(33) 

For tensor configurations 



rank 1 -» S'fr 1 0) = S[ Pft , S 1 (0\ fx) = S 2 p„, 

rank 2 -> S 1 (fi, v | 0) = S[ l2 5^ + S[ n p^ p v , S'(fJ, \ v) = S{ 22 8 MV + S' 121 p^p v , 

S J (0|/X,Z/) = S 222 8^, v + S221 PuPvi 



(34) 



with I = A,C, D, E, and we extract the scale dependence in complete analogy with Eq. 
Ultraviolet decompositions. The UV decomposition introduced in Eqs.© and (|7|) reads as 



S j( S ,{m}) - «j («,{"»};*) F 2 ( S ), 



(35) 



fe=-2 



where / = A, C, D, E classifies the family of integrals and j contains the information about the tensor 
structure. It is worth noting that the following results can be written in a compact form through the 
one-loop, two-point functions introduced in Subsection 13.21 We shall be quoting the tensor configurations 
relevant for our calculation. 



A family 



C family 



Sq (s, {m} 



123 : 



-2) = 2 



3 2 



S ( S > W 



123 • 



3 2 

-a s 



(36) 
(37) 



so 0> W 



1234 : 



Si (S, {m}l234 ; -2) 



(s,{to} 



1234 ■ 



-2) = 1, 



(38) 



s o (s,{™}i234; -1) = -26o(l,l, s,{m} 3i ; 0) - 1, 
s i ( s > {™}i234 ;-!)=- 6i(l, 1, s, {m} 34 ; 0) + -, 



s 2 ( s ,{ TO }l234 5 -1) 



2& 1 (1,1,s,{to} 3 4; 0) + 



1 



(39) 



D family 



s (s,{™}l2345 ; 


-2) 


= 0, 


1 

2' 




s 112 ( s ,{™}l2345 5 


-2) 




s 122 (s,{m 


s tjl ( 5 ,{™}l2345 5 


-2) 


= 0, 







-2) = 0, 



(40) 



8 



d Q ( S > {^}l2345 , -4.; — u, 

112 ( s > {™}l2345 ; -1) = - ^ b (l, 1, S, {m} 4 5 ; 0) 

1 



s 222 ( s > {™}l2345 5 ~1) 
s yi 0,{ m }l2345 5 -1) = 



1) = -i&o(l,M,{m}i2;0) 



3 

8 : 
3 



i? family 

S («:{™}l2343 
s 112 («,{™}l2343 ; -2 



2) = 0, sf (s, {m}i2343 ; - 

°> s 122 (S, {™}l2343 ; 



S ijl ( s , {™}l2343 i -2) = 0, 



-2) = 0, 



sf (s, {m}i2345 ; -1) = 0, 

s ?22 i s i {™}l2345 ;-!) = - -, 



1 

4' 



s 222 



(S,{m}l2343 ■ 



(41) 



1 

2" 



(42) 



s O ( s ; { m }l2343 ! 


-1) 


sf (s, {m}i2343 ; 


-1) 


s 2 ( s > {"l}l2343 ; 


-1) 


s 112 (s,{"l}l2343; 


-1) 


Sf u (s, {m}i2343 ; 


-1) 


s 122 ( s ,{ m }l2343; 


-1) 


s 121 ( s , {™}l2343 ; 


-1) 


s 222 ( s i {"2}l2343 \ 


-1) 


s 221 (S,{ m }l2343; 


-1) 



-2 6 (2,l,s,{m}34; 0), 
-6i(2,l, a, {m}a4 5 0), 
-2 6i(2,l,s,{m} 34 ; 0), 
1 m\ + m\ 
2 



6 (2, 1, s, {m} 34 I 0) + - 621(2, 1, s, {m} 34 ; 0), 
s 6 



-- 621(2, l,s,{m} 34 ; 0), 
-622(2, l,s,{m} 34 ; 0) - 
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-621(2, 1, s,{m} 3 4; 0), 
-2 6 22 (2,l,s,{m} 34 ; 0) - -, 

- 2 6 2 i(2,l, S ,{m} 34 ; 0). (43) 
The results shown for the s B family hold when m 3 7^ 0. 

Finite components. The UV- finite parts of two-point two- loop functions can be found in Ref. [10] . Eq. (89) 
and Eqs. (146-147), for S£ = 5 nl ; in Ref. [H], Subsection 5.8 for Sg = 5 121 , Subsection 7.3 for Sg = 5 131 
and Subsection 7.8 for Sfi = S$ 21 . 

Special configurations. Some conhgurations where one or more masses vanish show overlapping UV and 
IR divergencies and deserve a separate discussion. 

As a non trivial example, we consider the scalar function of the C family evaluated at p 2 = (sq ), since 
here the residue of the single UV pole can be IR divergent. Note that this function represents a tadpole, 
and it can be treated by means of the methods of Subsection 13.31 The residue of the double UV pole of Sq 
can be found in Eq. (|38]) and the residue of the single UV pole is given in Eq.(|39|). The same UV components 
arise when p 2 — and m 3 = m 4 is not vanishing, 



2 

s£(0, mi, m 2 ,m 3 , m 3 ) - -2F 2 _ 2 {M 2 ) - F 2 {M 2 )(l - 21n-g|) + s£ n , 



(44) 



where Sg n is finite when e — > 0. However, reduction in sub-loops requires to evaluate cases where m 3 = m 4 
0. Here the residue of the single UV pole becomes IR singular, and the explicit expressions read as 



s£(0, mi, m 2 , 0, 0) = 2 F 2 2 (M 2 ) + 2 

1 IR — I ^ <-*UV 

e \e 



( 2 , m l 2 1 m 2 \ 



F^(M 2 )+F 1R , 



(45) 
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Overlapping UV and IR divergencies appear also for special configurations of the scalar function of the E 
family (sf), evaluated at arbitrary s, when 7773 = 0, i.e. 



^(s,0,m 2 ,0,M,0) = 



l-x 



A 



-M-3 



i-A 

e 



2 ln(x - 1) + In 



Sq(s, mi, 7712,0, M, 0) 



1 -. 
In a; 

In x 



In a; 
_ ri 



- - A 



w] (t 1 

ln(ar-l)} + 



2 A t7V (s) - 1 - In a; 

S f n (s,0,m 2 ,0,Af,0), 
1 



2 hu>- 1) 
1b(x - 1)} - 



7715 



In- 



2A uv (s) - 1 



In- 



sf n (s,mi,m 2 ,0,M,0), 



(46) 



where Sg n is finite when e 



M 2 /s. Note that here double poles arise, whereas the residue of the 



and x 

double UV pole vanishes in Eq. ([42|) . 

These spurious singularities are a consequence of the sub-loop reduction of Ref. |19j . This technique 
allows to handle integrals with irreducible scalar products in the numerator and introduces new propagators 
with a corresponding zero mass. IR poles are then generated by the related integrals behaving as 

Jo 1 

The treatment of overlapping UV and IR poles is not free from ambiguities. Strictly speaking, the integrals 
of Eqs. pSj) and (|4"6"]) cannot be defined in any strip of the complex 77-plane. UV regulation requires n < 4 
and IR one should be performed for n > 4. Therefore, one should first of all disentangle overlapping UV and 
IR singularities. This can be achieved by means of the IBP identities introduced in Ref. [2D], which allow 
to write the original integral as a combination of objects which are UV or IR singular but never both. At 
the end one can use the well-known relation between UV and IR regulators, i.e. e = —e. We verified that 
this recipe gives the same result of the naive approach where one regularizes integrals for 77 < 4 without 
separating the UV and the IR domains. 

3.5 Two-loop three-point integrals 

Notation. Our notation for three-point two- loop integrals can be found in Section 9 of Ref. [14]. Here 
P = pi + P2 and we extract s = —P 2 in order to define scaled functions. For scalar configurations we 
introduce 



U 



-2 M.K.H 
Vn 



For tensor integrals, first we introduce form factors according to Section 9 of Ref. [T3]. Then, we define 
dimensionless functions in complete analogy with Eq. (|4"5)) . The (somehow redundant) definition of the 
rooting of momenta in two- loop vertices has been introduced for all families in Ref. [M]; we stick to the 
notation used there and one should only remember momentum conservation, i.e. P = pi + P2- 

Ultraviolet decompositions. The UV decomposition introduced in Eqs.© and Q reads as 



({momenta} , {m}) = Vj ({momenta} , {m}; k) F£(s), 



(49) 



fc=-2 



where L = E, I, G, M, K, H classifies the family of integrals under consideration and j summarizes the 
information about the tensor structure. Here we show only UV-divergent tensor configurations relevant for 
our calculation. 

E family 



V21 (• 



-2) = -2, 
-2) = 1, 



vfi (■ 

^2 (' 



-2) = 



2' 



J 12 1 



-2) = 1, 



-2) = 2, 



(50) 
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Vq (P2,P, {m}i2 34 ; -1) = -2&o(l, l,Pi,{m} 3 4; 0) - 1, 
ufi (P2,-P, {m}i234; -1) = -6i(l,l,pi,{m}34; 0) + -, 

v*2 (P2,-P,{m}i234; -1) =6o(l,l,Pi,{m} 34 ; 0) + 

«fi (P2,-P,{m}i234; -1) = -26i(l, l,pi,{m} 34 ; 0) + i 

«22 (P2,-P, {m}i23 4 ; -1) = 2&o(l,l,f>i,{m} 34 ; 0) + 



"114 (■• 



2) = 0, <■(■ 
2) = 0, 
-2) = 0, 



-2) = 0, 

-2) = 0, ; -2) = - 1 ^ 24 ("-;-2) = -^ 



w o(Pl> P >{ m }l2345; -1) = 
(pi,P, {m} 12345 ! -1) = 



C21 



w 2i (PIi-PjW} 12345 5 -1) = 
v[ 14 (Pl,P,{ m } 12345 5 -1) = 

(1,1,1, pi, P2,-P,{m} 345 ; 0) + 

+ 



^124 



-2c (l,l,l,pi,P2,-P,{m} 3 4 5 ; 0), 
-cu(l,l,l,pi,p2,P, {m} 34 5; 0), 
-2ch(1, 1, l,pi,p 2 ,-P, {m} 34 5; 0), 

1 THj^ 1 I 777 

- ^ 2 co(l,l,l,p 1 ,P,{m} 3 45;0) 

t C2 3 (l,l,l,Pi,P2,-P,{m} 3 

6 s 

1 p 2 

- — C22(l,l,l,Pi,P2,P,{m} 345 ; 0), 
2 

--C 2i (l,l,l,Pl,P2,-P, {"1)345; 0), 

-C24(l,l,l,Pi,P2,-P,{m} 345 ; 0) - 

-c 2 i(l,l,l,Pi,P2,-P,{m} 3 4 5 ; 0), 
„ , „ , , 1 



i (Pl,-P,{m} 12345 5 -1) = 

_ (pi,P,{m) 12345; -1) = 

^12* {Pl,P,{ m } 12345 5 -1) = 

w 224 (Pl,P,{m} 12345; -1) = -2c 2 4(l,l,l,Pl,P2,P,{m} 3 45; 0) - - 
^22i (Pl,-P, {»«} 12345 5 -1) = -2C 2 »(1, 1, l,Pl,P2,-P,{m} 3 45; 0). 



« G (• • • ; -2) = 0, «?• (• • • ; -2) - 0, 

«fu("-;-2) = -^ «?24 (' • • ; -2) = 0, ^ 24 (--- ;-2) = -\, 

v% k (• • • ; -2) = 0, 

<(■■■;-!) 

«n4(Pl>-Pl> P >{ m }l2345; -1) 
u ?2 4 (Pi i Pi ,P,{™} 12345; -1) 
«224 (Pl,Pl,P, {"l)l2 3 45 ! -1) 



-0, «5(...; -1) = 0, 

= - g b o(l,l,P2,{m} 45 ; 0) - -, 
1 

1 3 

= ~o b o(l,l,Pi,{m}i 2 ; 0)- -. 
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M family 

Vq (pi,P,{m} 123453; -1) = -2c (2, 1, l ) pi,p 2 ,P,{m} 3 4 5 ; 0), 
v ii(Pi, p A m } 123453; -1) = -cii(2,l,l,pi,p2,-P,{m} 3 45; 0), 
V2i(Pi, p ,{ m } 123453; -1) = -2ch(2, l,l,pi,p 2 ,-P, {m} 3 4 5 ; 0), 

1 YT~l^ j TYl^ 

V 1U (Pl, P ,{ m } 123453 ; -1) = X ~ -Co(2,l,l,pi,p 2 ,P,{m} 3 45; 0) 

2 s 

^ ^2 2^ P2 

+ c — c 2i(2, 1, l,Pi,P2,P, {m} 3 45 ; 0) + c 23 (2, 1, l,p 1 ,p 2 ,P, {m} 3i5 ; 0) 

b s b s 

1 p 2 

+ 7 — C22(2, 1, l,pi,p 2 ,-P, {™}345 ; 0), 

6 s 
2 

(Pl,-P,{»"}i23453; -1) = --^C 2 i(2,l,l,p 1 ,p 2 ,P,{m}u5; 0), 
2,:(Pi 5 ^{ m } 123453; -1) = -c 2l (2, l,l,p!,p 2 ,P, {m} 345 ; 0), 
2i(Pi>-P){ TO } 123453; -1) = -2c2i(2,l,l,p 1 ,p2,-P,{m} 3 4 5 ; 0). 



families 



v¥2i(Pi> p >{ m } 123453; -1) = -2c 2 j(2, l,l,p 1 ,p 2 ,-P,{m} 3 45; 0). (56) 
< (•••;-!) = 0, «5 (•••;-!) = 0, 

<14 (-P> Pi , P ,{ m } 123456! -1) = C o(M>l,Pl>P2,-P,{TO}456; )' 

vf jk {--- ;-l) = 0, '•;;,(•••: 1: (I. (57) 



4 M5 1 renormalization at two loops 

We now briefly summarize the notion of renormalization procedure. Any renormalization procedure for 
a gauge theory requires three steps: 

- the definition of the gauge-fixing term of the Lagrangian; 

- the renormalization-scheme prescription; 

- the choice of an IPS. 

In Subsection l4. II we define the gauge-fixing part of the Standard-Model Lagrangian and in Subsection l4.2l we 
discuss the renormalization scheme employed in this paper. This is a rather technical subject, but one should 
realize that choosing a renormalization scheme is equivalent to changing the variables which describe the 
degrees of freedom of the theory. Bare quantities are traded for renormalized ones through multiplicative 
renormalization constants, and the latter are expanded in perturbation theory introducing counterterms. 
The prescription imposed for the counterterms defines the so-called renormalization scheme. 

Whatever the renormalization scheme employed, a theory becomes predictive when renormalized quanti- 
ties are expressed in terms of the chosen IPS. As long as the same data are employed, different renormalization 
schemes lead to the same theoretical predictions. Special care should be devoted to the choice of the POs 
used as input: two predictions for the same PO differ by an amount proportional to the missing higher-order 
corrections if they refer to different IPS. We will discuss this issue in details in III, where we will show 
that a consistent description of unstable particles at two loops can be achieved in the framework of complex 
poles [21]. 



4.1 Gauge fixing 

We fix the gauge for the electroweak sector of the Standard-Model Lagrangian introducing six gauge 
parameters & (i = A, Z, AZ, tp°, W, <p), 



S*EW 

L 9f ~ 



c+c~ - - [ 

L L 2 



(c A y + (c z y 
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SA t,Z W 

Here j4 m , Z M and are the fields for the photon and the Z and the W bosons and <p° and ip^ are the fields 
for the neutral and charged unphysical Higgs-Kibble scalars. For the QCD sector, we employ the usual 't 
Hooft-Feynman gauge. 

Note that gauge parameters are usually introduced in order to test predictions for physical (pseudo- 
)observables (hereafter POs); S'-matrix elements are gauge independent, and one can indeed check that 
gauge-parameter dependence cancels out. Here, instead, the gauge-fixing term which defines what we call 
the R (( gauge is closely related to our renormalization strategy. We will discuss this issue in Subsection 14.21 
where we will introduce renormalized gauge parameters. 

Dyson-resummed propagators in the gauge. In this Section we illustrate the method of Dyson 
resummation of higher-order corrections which plays a crucial role in any renormalization procedure, because 
mass renormalization and wave- function factors for external legs naturally arise in this context. In the 
gauge, in addition, transitions between the vector gauge bosons and the Higgs-Kibble scalars, as well as a 
mixing between the photon and the Z boson, take place at the Born level. Therefore, Dyson resummation 
is a mandatory step already at lowest order in perturbation theory. 

Here we show the Dyson-resummed propagators for the vector gauge bosons and the Higgs-Kibble scalars, 
including all the allowed transitions. Furthermore, we provide the explicit expressions for the propagators 
of the Faddeev- Popov ghost fields. 

The charged sector: The Dyson-resummed propagator for the W boson reads as 

-jr-WW 1 1 / . -t-ww \ , 

A ^ = Wfi p^n^ O + App p » Pv )' (59) 

where we introduced 



-t-ww 1 



Moreover, for the propagator of the ip boson we get 



(P 2 + £w£ V M2)2 



(60) 



AW i p 2 + e w m 2 



Finally, a W-tp transition takes place, 



The neutral sector: Let us start with the Dyson-resummed propagators for the photon and the Z boson, 

= (2^i V 2 ( ^ + ^ v " Vv ) ' = WFi WTm^ ( ^ + s » p ^ ) ' (63) 



where we defined 



jZZ 

pp — 



V 2 + (p 2 + ^VM 2 f : 



+ & (p 2 + m 2 ) °" z r; ,/r:.L yQ - . (64) 



and we introduced a combination of gauge parameters, = £|/(1 + £| Clz)- Furthermore, the propagator 
for the <p° boson reads as 



1 ,r ■ (i. Ml 



(2n)H (p2 + 6L.M 2 f 
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In the neutral sector we have vector-scalar transitions, 

" (2^) 4 * (p 2 + 6 C v o M 2 ) 2 ^ " (27r)4t (p 2 + t z & M 2 ) 2 l&bj 

as well as a mixing between the photon and the Z boson, 

^ AZ i Lue z (p 2 + ^qM 2 ) 

~ (27r)«t p 2 (p 2 +£ z ^Af 2 )2 ^ 1 j 

Ghost fields: Finally, the gauge-fixing Lagrangian defines the propagators for the charged and neutral ghost 
fields X ± 1 Y z and Y A , 



~ (2tt) 4 z p 2 + £ v £w M 2 ' (2n)H p 2 + £ z M 2 ' (27r)*i p2 



4.2 Definition of the renormalization scheme: A/5 and beyond 

The orthodox approach to renormalization [5] uses the language of counterterms. It is worth noting that 
this is not a mandatory step, since one could write directly renormalization equations that connect the bare 
parameters of the Lagrangian to an IPS, skipping the introduction of intermediate renormalized quantities 
and avoiding any unnecessary reference to a given renormalization scheme. In this approach, carried on at 
one loop in [22], no special attention is paid to individual Green functions, and one is mainly concerned with: 

- UV finiteness of S-matrix elements after the proper treatment of external legs in amputated Green 
functions, which greatly reduces the complexity of the calculation; 

- gauge independence of predictions for (pseudo-)observables. 

However, renormalization equations are usually organized through different building blocks, where gauge- 
boson self-energies embed process-independent (universal) higher-order corrections and play a privileged 
role. Therefore, their structure has to be carefully analyzed, and the language of counterterms allows to 
disentangle UV overlapping divergencies which show up at two loops. In a renormalizable gauge theory, 
in fact, the UV poles of any Green function can be removed order- by-order in perturbation theory In 
addition, the imaginary part of a Green function at a given order is fixed, through unitarity constraints, 
by the previous orders. Therefore, UV-subtraction terms have to be at most polynomials in the external 
momenta (in the following, local subtraction terms). In this paper, we will express our results using the 
language of counterterms: 

- we promote bare quantities (parameters and fields) to renormalized ones; 

- we fix the counterterms at one loop in order to remove the UV poles from all one-loop Green functions; 

- at this stage, we check that two-loop Green functions develop local UV residues; 

- our final task is to fix the counterterms at two loops and to remove the UV poles from two-loop Green 
functions. 

Obviously, the absorption of UV divergencies into local counterterms does not exhaust the renormalization 
procedure, because we have still to connect renormalized quantities to POs, thus making the theory predic- 
tive. This will be the subject of III. In the remainder of this section we discuss renormalization constants 
for all parameters and fields. 
Renormalization constants 

We relate bare quantities to renormalized ones introducing multiplicative renormalization constants Zi and 
(if not otherwise stated) we expand them through the renormalized SU(2) coupling constant g R , 

oo / 2 \ n 

= * + £ if^J ^ (69) 
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where 8Z\ n are counterterms and the subscript i refers to masses, couplings, gauge parameters and fields. 
Let us make more precise the notion of renormalization constants. 

- for the bare masses of the W boson (M), the Higgs boson (M H ) and the fermions (to/) we write 

1/2 

m = Zm m Rl where m — M, M H ,m,f and m R corresponds to one of the associated renormalized 
masses. Z m is expanded through Eq. (frJ9")l . 

- For the bare SU(2) coupling constant g and the bare cosine (sine) of the weak-mixing angle eg (sg) 
we define p — Z p p R , where p = g, eg, sg and p R is one of the related renormalized parameters. Z p is 
expanded through Eq. (p)9"|) . 

- For the bare gauge parameters introduced in Eq. ([58|) we use £ = Zi=£ a , where £ is one of the bare 
gauge parameters and £ R is the associated renormalized quantity. Z% is expanded by means of Eq. (|69p 
except for the case £ = £ AZ , where we use 

n— 1 ^ ' 

After the expansions, we use the freedom in choosing the values for the renormalized gauge parameters 
and we set £ = 1, thus compromising between two alternative exigencies. On the one hand, countert- 
erms for the gauge parameters are required to remove UV divergencies; on the other hand, we recover 
part of the simplicity of 't Hooft-Feynman gauge. 

1/2 

- For a given bare field 4> it is convenient to write tfi — Z^ <p R , where 4> R is a renormalized field, and we 
expand Z^ through Eq. (j69|) . The bare photon field represents an exception, and here we use 




where A R and Z R are the renormalized fields for the photon and the Z boson. Note that Z AA is 
expanded through Eq. fM]) . 

In addition, bare fermion fields ip (we omit flavour labels) are written by means of bare left-handed 
and right-handed chiral fields ip L and ip R . The latter are traded for renormalized fields ip R and ip R 
expanding the renormalization constants through Eq. (|69[) . 

= \ (1 ± 7 5 )^, i> L ' R = ZT L , R ( ?2 ) 
Finally, Faddeev-Popov ghost fields are not renormalized. 



UV decompositions for Green functions 

Before defining our renormalization scheme, we introduce UV decompositions also for Green functions. 
Given a one- (i = 1) or two-loop (i — 2) Green function with TV external lines carrying Lorentz indices /ij, 
j = 1, . . . , N, we introduce form factors, 

A 

GL a„=VGliC uw , i = 1,2. (73) 

a=l 

Here the set K a , with a = 1, . . . , A, contains independent tensor structures made up of external momenta, 
Kronecker-delta functions, elements of the Clifford algebra and Levi-Civita tensors. Next, we introduce UV 
decompositions for the form factors G\ in analogy with Eqs.(H])-([7]), 

1 o 
1 loop - G\ = Gl ]k F^M 2 R ), 2 loops -> G 2 a = £ G 2 a . k F 2 (M R ). (74) 

k=-l k=-2 
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The UV factors (k = —1,0, 1) and F% (k = -2,-1, 0) can be read in Eqs.® and ©. 
Renormalization prescriptions: A/5 and beyond 

In the spirit of the UV decomposition of Eq. ([74|) . one could define a non- minimal (NMS) subtraction scheme 
where the one- and two-loop counterterms are defined by 

-l 

1 loop - = AzPf^MD, 2 loops -> 5z\ 2) = ]T Az\%F 2 {Ml), (75) 

fe=-2 

and are fixed in order to remove order- by-order the poles at e = for any Green function. As a result, also 
the product of a one- loop counterterm with a one- loop diagram (i.e., a one- loop counterterm insertion) has 
the same UV decomposition of a two-loop function thus simplifying two-loop renormalized Green functions. 
Even if we scale non-tadpole diagrams with s instead of M|, the difference is proportional to In M|/(s e), a 
non-local term which cancels in the total. 

Although the NMS scheme has the virtue of respecting a universal UV decomposition, we adopt the 
conventional M S choice, keeping the counterterms as simple as possible. Let us accordingly define 

SZP = (-7 + c) A ^ (1) ; (76) 

and we define the renormalization scheme choosing the explicit value for C. We define a minimal MS 
subtraction scheme, 

SZ^ = f-^+AuJ) Azf, (77) 

and we fix the counterterms in order to remove the poles at e = for any one-loop Green function. Next, 
we extend the scheme at two loops by introducing 

«tf°= (^ + AZg?) +Al v AZ$. (78) 

Accordingly, we decompose the two-loop form factors, 

G 2 a = — J G 2 a . uv . x + (- — Auyj Gl. uv . 2 + A| v G 2 a . uv ;3 + G 2 ■ p , (79) 

and we fix the counterterms in order to remove the poles at e = for any two-loop Green function. Note 
the correspondence with the decomposition in Eq. l|71|) , 

W ; UV ; 1 ~ W ; -2i lj a;[Jl';2 _ "a;-li 

G 2 a ■ uv ;3 — 2 I -2' G 2 a . F — G 2 . . (80) 

5 One-loop renormalization 

In this section we collect all the one-loop results necessary for our renormalization procedure. Tadpole 
renormalization (/? vertices) and neutral-sector diagonalization (r vertices) are shown in Subsection 1 5 . 1 1 and 
counterterms for parameters and fields are presented in Subsection 15.21 

In Subsection 15 . 31 we discuss the relevance of one-loop finite renormalization through some examples and 
in Subsection 15 .41 we analyze how counterterms affect WST identities. Renormalized WST identities riaffcrm 
the basic simplicity of the one-loop structure of the theory. These results should provide a useful guide to 
the two-loop renormalization. 
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5.1 P and T at one loop 

Two key ingredients are the parameter /3 t , introduced in Sect. 2.3 of I (tadpole renormalization), and the 
parameter T, defined in Sect. 3.1 of I (all-order diagonalization of the neutral sector). We recall that one 
can define different tadpole schemes, and in I, we discussed the so-called flh and (3 t schemes; since here no 
ambiguities arise, we drop everywhere the subscript t from (3. To show the one-loop results, we expand /? 
and r in the unrenormalizcd SU(2) coupling constant g, 

oo oo 

= Y,9 2n Pn, r = ]T 3 2 "r„. (8i) 

n— 1 n—1 

The complete answer for (3\ reads as 

^3 3 

Pi = 16n 2 x {^^2[ x d,i a o( m d,i) + a: «,» Q o( m »,')] + ^2xl t a (mi^) 

H 2—1 i— 1 

- ^1 (^T 1 + ^)ao(Mo) - \(n- 1 + ^)a (M) - ^4«o(Mj}, (82) 

9 

where we recall that Xi = mf/M 2 , li is the charged lepton and Ui, d\ are the up and down quarks of the ith 
fermion doublet. The one-loop tadpole ao is given in Eq. (fT3|) and the explicit result for Ti is 

n — 2 

ri= 16^ ao(M) - (83) 

Note that in Eqs. (f8"2"|) and ([55)1 we dropped the subscript R, since bare quantities disappear after the intro- 
duction of renormalized parameters and fields. 

5.2 One-loop counterterms 

In this subsection we provide the full list of one-loop counterterms in the MS scheme. Counterterms 
arc defined in Eqs. (j6T))) - (jT2")) , and UV factors are extracted through Eq. ([77)) . Here we introduce short-hand 
notations for sums over fermions (I — > charged leptons, u,d — > quarks), 

i—l i—1 2—1 

Gauge parameters. Gauge-parameter counterterms can be expressed by means of field, mass and coupling- 



constant counterterms, 
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- 1 A7 (1) 


AZ^ = 


40 



(1) 

AZ ? 



4 +AZ«j. (85) 

The gauge-fixing term defined in Eq. (|58p is invariant under renormalization. We will discuss this issue in 
Subsection [ 



Gauge-boson and Higghs-Kibble fields. 



^21 = 2 4sl AZ^ = -^--23ce), 



3 3 V 2 c e 

1 /41 1 . \ . „m . 1 



AZ ^ = L ( li -L _ 41 + 23 4) , AZ« = -1 - i 

3V2 c V ' v° 2 



-o — Xi — 3 ( + X u 



AZ$ = \, AZW=AZ«- (86) 
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Masses and couplings. Here we discover a nice feature of the neutral-sector diagonalization introduced in 



Sect. 3.1 of I, since in the full-fledged Standard Model we derive a QED-like relation involving AZ^\ AZ£ 
and AZ { 2}, 



AZ 



(i) 



3 1 
-14 
1 /41 1 
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3 



3 1 

2cl 



AZc » -2U d 2 + 3 C< 



3-2X ; 2 -6(x d 2 +^) 



-re* + X, + 3 ( X d 



X 



AZ^ = -AZ^- 1 -AZ { }1 (87) 



In other words, electric-charge renormalization depends only on the photon vacuum-polarization function, 
also in the (complete) Standard Model. 

Higgs-boson field and mass. 
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A7 (i) _ 6 



11 1 1 v 

T7 + 1 X„ Xl 

2 c 2*3* 
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Fermion fields and masses. Here v is the neutrino field associated with the charged lepton I, and u and d 
are up and down quarks belonging to the same doublet. Furthermore, g s is the SU(3) coupling constant of 
strong interactions. For left-handed spinor-field countcrtcrms we get 



a m 1 Z 1 1 1 
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(89) 



Expressions for the right-handed components are given by 
AZ" - 0, 





4 


1 


4 


«s 


~ 9 


r 2 


~ 9 


AZ« 


1 


1 


1 


OH 


~ 9 


r 2 


~ 9 



2 " 
1 

2*" 



4ff| 
3 5 2 ' 

3 ff 2 ' 



(90) 



Finally, fermion-mass counterterms depend on the sum over the fermion families because of tadpole renor- 
malization (the explicit expression for (3\ can be found in Eq. (|82p ). 
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(91) 



5.3 Finite renormalization 

The last step in one-loop renormalization is the connection between renormalized quantities and POs. 
Since all quantities at this stage are UV-free, we yerm it finite renormalization. Note that the absorption 
of UV divergencies into local counterterms is, to some extent, a trivial step; finite renormalization, instead, 
requires more attention. For example, beyond one loop one cannot use on-shell masses but only complex 
poles for all unstable particles. The complete formulation of finite renormalization will be given in III. 
However, let us show some examples where the concept of an on-shell can be employed. Suppose that we 
renormalize a physical (pseudo-)observable F, 



F = F B + g 2 F lL (M 2 ) + g 4 F 2l (M 2 ), 



(92) 
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where M is some renormalized mass which appears at one and two loops in F\ L and F% L but does not show 
up in the Born term F B . In this case we can use the concept of an on-shell mass identifying M = M os for 
the two-loop term and performing a finite mass renormalization at one loop, 



M 2 =M 2 oa \l 



9 



16 7T 2 



ReS 



SZ 



(i) 



}=M 2 s +g 2 AM 2 



(93) 



where M os is the on-shell mass and £ is extracted from the required one-particle irreducible Green function. 
Eq. (j9"3")) is still meaningful (no dependence on gauge parameters) and will be used inside the one- loop result, 



F = F B + g 2 F lL (M 2 s 



where 
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F 2L (M 2 S ) + F{ L (M 2 S ) AM 



dF lL (M 2 



g M 2 \M*=Ml s - 



(94) 
(95) 



Here we show some examples of one-loop finite renormalization. 
W boson. For the W boson we write 

,-,2 



M 2 = Ml 
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16tt 2 



Re Y> ww — SZ 



(i) 



(96) 



where the expression for the counterterm can be read in Eq. (|87p . The quantity within square brackets in 
Eq. (|96p is finite by construction, and it is convenient to write the decomposition 



E w „, + 



2(/3 1 + r 1 ). 



(97) 



Pi and Ti can be read in Eq. ([82|) and Eq. (|83| . whereas the first two terms in the r.h.s. of Eq. ([97| refer to 
fermionic and bosonic contributions. The fermionic part for one generation (one leptonic family, xi = mf / 
M 2 , and one up-down quark doublet, x u = rn^/M 2 , Xd = m 2 d / M 2 ) is 
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(98) 



and the bosonic component reads as follows, 
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Fermions. For a fermion / we obtain 
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(100) 
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where the fermion-mass counterterms can be read in Eq. (|91[) . and the physical mass in Eq. (|100p coincides 
with the on-shell one. Using 



Qi — — 1) Qu — g : Qd 



Nf = 1,^ = 3, v f = lf ] -2Q f s 



(101) 



(3) 

where Ij is the third component of isospin, we provide the needed expressions for arbitrary up and down 
partners in a given fermion doublet. For a generic up- type fermion we get 
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(102) 



For a d-type fermion, the result reads as 
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Higgs boson. We also provide the result for the Higgs-boson mass, 
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(104) 



where the counterterm can be found in Eq. (|88p . and we use a decomposition similar to Eq. (|97p . T, HH = 
T,{j H + S^ H . For one generation of fermions we get 
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(105) 



f—l,u,d 11 f=l,u,d 

whereas the bosonic component reads as follows, 
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b (l,l,M 2 H , M,M) 
- (-^r + -yjj 6o(l, 1 , M 2 , M , M ) + -x%b (l, 1,M 2 H1 M H ,M H ). 



(106) 



The sum over the fermion masses in Eq. (|106p depends on the (3\ parameter (see Eq.(|82p). which we included 
explicitly in the computation of the Higgs-boson self-energy. 
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5.4 Ward- Slavnov- Taylor identities 

Ward-Slavnov- Taylor (WST) identities play a crucial role in showing that the electroweak theory is 
renormalizable and respects unitarity at the same time. In fact, we can prove with their help that transition 
amplitudes are the same in every non-singular gauge. Having derived the one-loop counterterms, we analyze 
how renormalization affects the WST identities of the theory; note that here we consider only a particular 
class of WST identities, the so-called doubly- contracted ones with two external gauge-boson fields. Let us 
consider the gauge-fixing functions C a of Eq. (|58j) written in terms of bare quantities, 

C A = - t z d,Zl C z = ~d»Zl + ^oM V>°, C± = —d^W^ + $My ,:t . (107) 

SA Sz Sw 

Then, we introduce renormalized quantities and we set the renormalized gauge parameters to unity, 

C A = —ZXaZ^^B^A^ — (zJ z Z^ + Z£ Z£ Az ^d f _ l Z^, 

c z = -zlz-^d^z^ + z^zhz-iz^Mov , 

C± = -Z&Z^dpWjt + zlzlZ^M^. (108) 

Next, we associate a set of source vertices to the Fourier transform of a gauge-fixing function C a . 

Let us consider C , where we get the two source vertices of Fig. [T] A doubly-contracted two-point WST 
identity is obtained by connecting two sources through vertices and propagators. Here we get, at every order 
in perturbation theory, the identity of Fig. [2j We analyze the WST identity at O (g 2 ) . After expanding 

w 

p M 

# . . . (2^)4, MZ^Z^Z^ 

V 



Figure 1: Sources related to the gauge-fixing functions defined in Eq. (|108[> . The momentum p is flowing inwards. 



W W W 



W V _^ 



■ • =0 



Figure 2: Doubly-contracted WST identity with two external C sources. Gray circles contain all the irreducible 
and reducible Feynman diagrams contributing to the needed Green functions. Black dotted circles represent sources. 
Their expression can be read in Fig. [T] 



the renormalization constants Zi up to O (<? 2 ) we obtain 



-^-J-sztf 

167r 2 V2 



SZ. 



+ Uz^ ] +5Z^Mcp ± . (109) 



and we get the source terms of Fig. [3] As a consequence, the doubly-contracted WST identity with two 
external C ± sources at O (g 2 ) of Fig. [2] receives two different contributions, whose graphical representations 
are given in Fig. 3] In the first line we have Green functions containing Feynman diagrams at O [g 2 ) and 
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source vertices at O (1). The second line, instead, contains Green functions with Feynman diagrams at O (1) 
and source vertices at O (g 2 )- Finally we use Eq.i 



sz 



(i! 



(110) 



and we are immediately left with the four contributions displayed in the first line of Fig. |U Note that 
this is exactly the doubly-contracted identity with two external C ± -source vertices which appears in the 't 
Hooft-Feynman gauge. In order to get the complete set of doubly-contracted WST identities we insert the 
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Figure 3: Source vertices related to the gauge-fixing functions C . We start from Fig. [Tj and expand the renormal- 
ization constants up to O (<? 2 )- The black squares denote source vertices at O (1), whereas the white squares denote 
source vertices at O (g 2 )- 
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Figure 4: Doubly-contracted WST identity with two external C sources at O (g 2 ) . Gray circles contain all the 
irreducible and reducible Feynman diagrams contributing to the Green functions at O (<? 2 ) . Black squares and white 
squares represent source vertices. Their expression can be read in Fig. [3] 

explicit expressions for the one-loop countcrtcrms of Eq. ((85|) in Eq. (|108|) . Besides Eq. (|110j) we get 



bzf } = -6Z^>, SZ' 1 ' = -5Z)ti, 6Z' X > = -SZ", 1 ', 5Zf> = 



,(1) 17W 



2=^)-^+^). (in) 

Using Eq. lfTTUj) and Eq. (fTTTj) in Eq.([l08]) we finally obtain 

C A = -d^ + O {g A ) , C z = -d^Z^ + M Q ^ + 0{g A ), C± = -d^W* + Mcp^ + O ( 5 4 ) . (112) 

The whole set of one-loop WST identities is given in Fig.fSJ whereas source terms can be read in Fig. [5] The 
results are exactly the doubly-contracted identities with two external source vertices in 't Hooft-Feynman 
gauge. 
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Figure 5: Doubly-contracted WST identities with two external gauge bosons at O (g 2 )- Gray circles denote the sum 
of the needed Feynman diagrams at O (<? 2 )- Source vertices are given in Fig. [II] 
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Figure 6: Renormalized source vertices for the WST identities of Fig. [5] 



6 [3 and T at two loops 

The first step in extending our renormalization procedure beyond one loop has to do, once again, with 
tadpole renormalization ((3, Subsections 16. II and I6.2p and neutral-sector diagonalization (T, Subsections I6.3I 
and E3J). 



6.1 fc: reducible contributions 



We consider the contributions of Fig. [7J to the reducible Green function with one external Higgs boson 
atO(.g 3 ), 
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(113) 



where the sum runs over the nine families of Feynman diagrams of Fig. [7] The sum over the first six families 
is given by 
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(114) 



where F".^. denotes one of the three families of Feynman diagrams which contribute to the Green function 
with two external Higgs bosons at O (g 2 ) and F^.' iTT is one of the two families of Feynman diagrams which 
contribute to the Green function with one external Higgs boson at O (g). As a consequence of our choice for 
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Figure 7: The reducible contributions to the renormalization of the Higgs-boson vacuum expectation value at O (<? 3 ) . 
Nine families of Feynman diagrams give a contribution. Dashed lines represent a Higgs boson. Dotted one-leg vertices 
are /3i-dependent vertices at O(g). Dotted two-leg vertices, instead, depend on 0i and on one-loop counterterms. 
Combinatorial factors for two-loop reducible diagrams are given by the products of those for one-loop diagrams, 
except for two cases, where we indicated explicitly the additional multiplicative factor. 



Pi we get 

i=i 

For the last three diagrams of Fig. [7] we obtain another vanishing component, 

EpH-2 _ 1 3 gM| 1 r a P ff;i] 2 _ n nir x 
7 ^;rcd" (2 ^4 i4 M ( p 2 +M 2)2 ^1 ;irr + ^2 ;irrj U ' 

As a result, Pi receives contributions only from irreducible diagrams. 

6.2 p2- irreducible contributions 

In order to derive P% we take in account only irreducible diagrams and we write 

i=i 

where the sum is over the six families of irreducible diagrams of Fig. [5] Note that here we are dealing 
with three kinds of diagrams. The first three families contain two-loop diagrams at O (g 3 ) . The fourth and 
the fifth family contain one-loop diagrams obtained with the insertion of a dotted vertex which represent 
a Pi, T\ or one-loop counterterm insertion. Explicit expressions for Pi and Fi can be read respectively in 
Eqs.([52]) and |53]). One- loop counterterms are given in Eas.((55 )l -((55 ]) . Eqs.([S!l])-([II]) and Eq.®. The last 
family contains one diagram with a /32-dependent vertex. Here we present the UV divergent residue for Pi 
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Figure 8: The irreducible contributions to the renormalization of the Higgs-boson vacuum expectation value at 
O (g 3 )- Six families of Feynman diagrams give a contribution. Dotted vertices depend on /3i, Fi or the one-loop 
counterterms. The black-square vertex, instead, represents a /32-dependent vertex. 



neglecting the fermion masses, except for the top-quark mass m t . Although not reported in this paper, the 
full result is available. As usual, we write the decomposition 
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The expression for 02;uv,i is given by 
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The expression for P2;uv,2 reads as 
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Now, the expression for 02 ■ uv,3 is given by 
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It is worth nothing that, as a consequence of our choice for 02, reducible Feynman diagrams containing a 
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Figure 9: O (<? 3 ) tadpole insertions contributing to the reducible two-point Green function at O (g k )- The gray circle 
denotes the three-point Green function, with an additional Higgs boson line, at O (g k ~ 3 )- These diagrams cancel as 
a consequence of our choice for fe. Dotted and squared vertices are the same introduced in Fig. [8] 

tadpole subdiagram at O (g) 3 should not be included in any computation (this is displayed in Fig. [5]). 

6.3 T2: reducible contributions 

Let us start with the reducible contributions of Fig. 1101 whose contribution to the Green function with one 
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Figure fO: Reducible contributions to the Green function with one external photon and one external Z boson at 
O (g 4 )- Gray circles denote the sum over the irreducible Feynman diagrams at O (g 2 )- 

external photon and one external Z boson at O (g 4 ) can be written through the two-point Green functions 
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at one loop involving neutral gauge bosons, 



G 



Az ; 2 



[a> ; red " 

We extract form factors as 



| \ (~lAA\l 1 (~,AZ;l .f~,AZ\\ ' f^iZZ : 1 .s-yAip°;l j- ^! / 3°Z;l\ 

7r) 4 i I MCt ; '"p 2 ; irr Ma 1 ilr p 2 + Mq ; irr M ; irr p 2 + Mq ' irr J 



= Gj' 2 {p 2 )5 l _ iv + G^^PuPu vectors), 
= G l J ,2 (p 2 )ip^ (i vector, j scalar) 



and obtain 



M 2 



The photon self-energy shows factorization, 
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It follows that Eq. (|124p . at p 2 = 0, can be written as 
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Following our choice for Ti we haveG^'^O) = 0. As a consequence, the reducible contributions of Eq. (|124[) 
vanish. Therefore, T2 is fixed through irreducible Feynman diagrams, 
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6.4 r 2 : irreducible contributions 

To fix T2 we use 
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where the sum is over the form factors of the thirteen families of irreducible diagrams of Fig. 1141 Here we 
get three kinds of diagrams. The first eight families contain two-loop diagrams at O (g 4 )- Next, we have 
four families of one-loop diagrams obtained by inserting a vertex which depends on T\, or on a one- loop 
counterterm. Note that these diagrams depend also on the counterterms for the gauge parameters 5Z^ . 
Obviously, there will be no track of the gauge parameters in the final results since we made the identification 
£i = 1. Nevertheless, gauge-parameter counterterms play an essential role in making every Green function 
UV-finite. Finally, the last family contains a single Feynman diagram with a r2-dependent vertex. 
The contribution of the thirteenth family of Fig. Q3]is given by 
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Using Eq.(dl§ in Eq.{T2S|| we get 
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T 2 can be expressed through a linear combination of products of one- and two-loop scalar integrals, evaluated 
at zero-momentum transfer. As usual, we write the decomposition 



t 2 =[\-a uv ) ( T2 ' uvA + r 2;[7V , 2 ) + A 2 uv r 2 . uv , 3 + r 2 ., F , inn 



where the UV residues are 



169 137 111 



F2 ; UV ' X ~ 1536tt4 ' F2 ; uv - 2 ~ 6144tt4 ' r2 ; uv ' 3 ~ 2048*4 ' (132) 



7 Two-loop self-energies 

In this section we will show that two-loop self-energies develop UV poles whose residues are polynomials in 
the external momentum (local residues). Therefore, a suitable choice for the counterterms in the MS scheme 
removes UV divergencies also at two loops. We generate the needed Feynman diagrams with QraphShot [8] 
and express every two-loop two-point Green function through scalar integrals using standard tensor-reduction 
techniques. Finally, we extract UV residues by means of the results collected in Section [2l 

In Subsection 17.11 we introduce our notations, and in Subsection 17.21 and Subsection 17.31 we review the 
transverse form factors for the photon self-energy and the transition between the photon and the Z boson. 
Explicit results for the transverse components of the Z- and VF-boson self-energies and for the Higgs-boson 
self-energy are shown in Subsection 17.41 Subsection 17.51 and Subsection 17.61 

7.1 Notation 

Reducible two-loop two-point Green functions with two external boson lines are displayed in Fig. 111! 
They are UV finite by construction, since one-loop counterterms have already been fixed. Therefore, we can 
limit our discussion to irreducible two-loop two-point Green functions, which can be written through the 
form factors introduced in Eq. (|123j) . 

^;\l(s) = J2 F a%-,k(^ (133) 

fc=l 

Here i and j are the external fields, the subscript a labels the form factor and the sum runs over the thirteen 
families of Feynman diagrams of Fig. [TU We shall organize the list of results according to the following 
decomposition: 

^a;'irr( s ) = ^a;'irr ° S ( S ) + ^a-itv ^ ' ( s ) + 'j'irr ^( s ) + ^aj'irr * ( s )) (134) 

where we split the complete result into purely-bosonic diagrams (bos), diagrams with one lepton sub- loop 
(lep), diagrams with one light-quark sub-loop (lq; u, d, c and s) and diagrams with one top-bottom sub-loop 
(tb). Furthermore, we extract UV poles using Eq. (|T9")) and we always leave out an overall factor (27r) 4 i. In 
the following we show results neglecting the masses of the leptons and the five light quarks. Note, however, 
that the full result is available. 



7.2 The photon self-energy 

When computing the photon self-energy at a given order in perturbation theory we adopt the convention 
of excluding diagrams which contain light quarks coupled to photons or gluons with no other hard scale. 
In this case the correct treatment requires the use of non-perturbative dispersion relations [23 . So there 
are two preliminary steps: first, we exclude from the computation two-loop diagrams which contain only 
light quarks and photons (or gluons). Next, we derive again one-loop counterterms for light-quark fields and 
masses neglecting one- loop diagrams which contain only light quarks and photons (or gluons). The results 



28 



Q Q 



n 



n 



n 



o- 



Figure 11: Reducible contributions to the two-point Green function with external boson fields i and j at O (p 4 ) 
The dotted two-leg vertex represents a /3i-, IV or one-loop-counterterm-dependent vertex. 



for up- and down-quark fields reads as 
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Light-quark mass counterterms are given by 
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Two-loop result 



Before showing explicit expressions for the individual components of Eq. Q134p let us recall that scaled masses 
arc defined by Xi = M 2 /M 2 , where M is the renormalized M^-boson mass, and where yt = Mf/s, with 
s = —p 2 ; p is the external momentum. Each term in the r.h.s. of Eq. (|134p develops a non-local residue 
containing the (logarithmic) function 

L (x) = In §4^' ( 137 ) 

(3(x) - 1 

where f3 can be found in Eq. (|21[) . We will show explicitly that non-local residues cancel when summing up 
all the contributions, and UV poles will be subtracted by polynomial counterterms. 

Two light- quark doublets 
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Top-bottom doublet 
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Note the presence of mixed electroweak-QCD corrections proportional to the strong coupling constant g s , 
which do not show up in Eq. (|139j) because they are subtracted from the light-quark components. 



Three lepton doublets 
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Bosonic contributions 
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Total contributions 



The sum of the four components gives 
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Non-local residues cancel and the UV residue can be reabsorbed by the counterterm contributions, 

Cd;irr;'[/v;a( S ) = oi:fi_4 i «' ° = I' 2 ' 3 - 
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An alternative method 
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Figure 12: Counterterm-dependent Feynman diagrams with a y> loop. Dotted vertices represent counterterm 
insertions. 

Here we outline another strategy to verify the cancellation of logarithmic residues. Let us consider Fig. [T2l 
corresponding to two Feynman diagrams with a Higgs-Kibble ^ loop and a one-loop counterterm insertion. 
The contributions of these two diagrams is given by 
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After summing Eqs. (|151[) — (|152p the terms which depend on 8Z^p cancel out, 
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Moreover, after using the relation 
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the answer contains only the VT-mass- and gauge-parameter-dependent counterterms, 
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Furthermore, all Feynman diagrams with the insertion of counterterms can be paired like in Fig. 1121 and we 
can repeat the argument for any of the pairs. Consequently, counterterms for masses and gauge parameters 
suffice in removing logarithmic residues. Here we outline the procedure: we take in account the irreducible 
Green function at O (<? 2 ) which are written as follows: 



GAA ; 1 ; fcr / \ . s~fAA ; 1 : bos / \ 
d;irr ( S ) + G d : far ( S )> 



(156) 



where the first term receives contributions from the three lepton families and the top quark (light-quark 
components are subtracted), whereas the second one gets contributions from bosonic loops. The essential 
point is that only bare quantities appear in Eq. (|156p : in other words, the gauge-fixing terms are 



c- 



Sz 



d M W±> b + &M b <p k > b , (157) 
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where we used the superscript b to denote bare quantities. Next, we use fermion mass renormalization, 
m b = Zm 2 m, where the renormalization constant is expanded as 

Z m = l + ^5Z$ + 0{g A ), (158) 

and the counterterms for lepton and top masses can be read in Eq. (|9ip . We rewrite the first term in the 
r.h.s. of Eg. ([156)1 as 

GS# !fcr « = GSf&^M + (159) 
where functions in in the r.h.s. of Eq. (|159p depend on renormalized quantities. The first term in r.h.s. of 
Eq. (|159p represents the one- loop fermionic contribution to the renormalized photon self-energy, whereas the 
second term gives an additional contribution to the two- loop result. 

We verified that, after adding the second term in the r.h.s. of Eq. (|159p to two-loop diagrams, the answer 
is free from logarithmic residues. Therefore, fermion-mass renormalization removes all logarithms in the 
residue of the simple UV pole for the fermionic part. 

However, a non-local residue remains in the bosonic part. Unfortunately, a simple procedure of IV-mass 
renormalization is not enough to get rid of logarithmic residues in the bosonic component and the reason is 
that in a bosonic loop we may have three different fields, the IV, the p and the charged ghost fields, and 
only one mass is available. The situation is illustrated in Fig. [13] where the black dot denotes insertion of a 
counterterm SZ M . The latter is fixed in order to remove the UV pole in the W self-energy and one verifies 
that the total in the second and third line of Fig. [T3] (ip and X self-energies, respectively) is not UV finite. 
The procedure has to be changed if we want to make the result in the bosonic sector as similar as possible 






Figure 13: W-mass counterterm insertion in the charged one-loop transitions. The W — W one is UV finite, whereas 
the same is not true for ip — ip and ghost-ghost transitions. 

to the one in the fermionic sector. A diagrammatic interpretation of the method is displayed in Fig. 1151 For 
the bosonic part we have to take into account Dyson-resummed propagators. Let us consider the following 
integral, corresponding to a ip loop in the photon self-energy, 

F - -<*)mi> j d « w^m^w^m] (2? " +p " )(29 " + ^- (160) 

We introduce renormalized quantities for the gauge parameters and the TV-boson mass, 

M b = zi /2 M, ^ W =Z^ W , i p = Z^i„, (161) 
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and expand the renormalization constants, 



Zi = l+ if^f > + O (g 2 ) . (162) 



For the propagators we get 



[q 2 + = (q 2 + M 2 )- 1 - ^(SZ £w + 6Z 6v + SZ M )M 2 (q 2 + M 2 y 2 + .... (163) 

Therefore, each one-loop bosonic diagram can be written as 

BftZi') = H^.^is) + A//.-;: .„», (164) 

where the second term in the r.h.s. of Eq. (ll64|) gives an additional contribution to the two-loop photon 
self-energy which depends on the VF-mass and gauge-parameter counterterms. Collecting all the diagrams, 
we get 

GAA ; 1 ; bos/ \ /~iAA ; 1 ; bos/ \ . a /~tAA ; 2 ; bos/ \ /-. ar\ 

jr. in, (S) - 0- d . irl . . ron (S) + ^. irr . rcn (S). [LOO) 

We verified that PU-mass and gauge-parameter renormalization removes all the non-local residues also for 
the bosonic component of the photon two-loop self-energy. 

We emphasize the importance of using the gauge of Eq. (|58|) : renormalization of gauge parameters 
is essential in removing UV divergencies at one loop; furthermore, unitarity requires that any subtraction 
term must be local. 

Summarizing, we have been able to verify that the electroweak theory can be made (two-loop) UV finite 
by adding local counterterms with two different methods, each leading to the same result. In other words, 
the well-known one-loop result that self-energies suffice in performing renormalization can be extended up 
to two loops. Although counterterms have been explicitly included it remains true that: 

- In extracting a from Thomson scattering at zero momentum transfer we find four classes of two-loop 
diagrams: 

I) irreducible two-loop vertices and wave-function factors, product of one-loop corrected vertices 
with one-loop wave-function factors; 

II) one-loop vacuum polarization ® one-loop vertices or one-loop wave-function factors; 

III) irreducible two-loop A A, AZ, A<fi° transitions; 

IV) reducible two-loop AA, AZ, A<fi° transitions. 

We have verified that the non-vanishing contribution originates from III and IV only and, within these 
terms, only the reducible and irreducible AA transition survives. For this result the role of T is vital. 

- In extracting the Fermi coupling constant from the muon lifetime all corrections to 

G F g 2 



s/2 8M 



(l + Ag) (166) 



which do not originate from the W self-energy and that are UV (and IR) finite at one loop remain 
finite at two loops after one- loop renormalization (i.e. two- loop counterterms are not needed), 

The proof has been obtained by using QraphShot generating the whole set of corrections. The result 
follows by algebraic methods (see Appendix [Dj , i.e. full reduction of tensor structures without using the 
explicit expressions for the scalar integrals. 
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7.3 The 7-Z transition 

For the transition between the photon and the Z boson we show again all the various components of 
Eq.HSl. 

Two light-quark doublets 
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Top-bottom doublet 
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Three lepton doublets 
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Bosonic contributions 
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The sum of the four components gives 
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and can be absorbed by the contribution of the counterterms, 
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7.4 The Z-boson self-energy 

The explicit expressions for the individual components of the i?-boson self-energy are rather lengthy. Here 
we provide the total result for the coefficients of the UV factors and we show that non-local residues cancel, 
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As evident, the two-loop counterterms contain only local residues 
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7.5 The VK-boson self-energy 

For the VF-boson self-energy we provide the total result which shows once again absence of non-local residues. 
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7.6 The Higgs-boson self-energy 

Here we provide the result for the Higgs-boson self-energy, showing the absence of non-local residues, 
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Therefore, we can employ a non-local subtraction term, 
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8 Two-Loop Counterterms 

In this Section we provide the full list of two-loop counterterms (for reason of space fermion masses other 
than the top-quark one are not shown, except for the first result). 

Let us start with field counterterms, where the sum over the fermion masses was defined in Eq. (|84[) . 
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For the cosine of the weak-mixing angle we obtain 
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Counterterms for the Higgs-boson mass read as 
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Here we present also the two-loop counterterms for the Higgs-Kibble tp and tp scalar bosons, 
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The equalities in Eq.flSS]) and in Eqs. (j2"T4"|) - (f2"T7jl) are related to the custodial symmetry SU(2) L ® SU(2) R -> 
SU(2) V . Finally gauge-parameter counterterms can be expressed by means of Eqs. (|192|) ~ (|216|) . similarly to 
the one-loop case, 
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After deriving the relations for the two-loop gauge-parameter counterterms, we can verify immediately an 
important result: the doubly-contracted WST identities with two external gauge-boson fields have the same 
form derived in 't Hooft-Feynman gauge. The gauge- fixing functions given in Eq. (|58p are expressed by 
means of bare quantities. After introducing renormalized quantities, the gauge-fixing functions change; their 
expression, order-by-order in perturbation theory, can be read in Eq. (|108[) . By expanding Eq. (|108p up to 
two loops, and using the expression of Eqs. (|217|) - (I218|) for the gauge-parameter counterterms, we derive 

C A = -d M Ap + O (g 6 ) , C z = -d fl Z„ + M $° + O(g 6 ), C ± = -0 M W M + M$ + O (/) 

and the doubly-contracted WST identities are formally equivalent to those of Fig. at O (g 4 ) . 
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9 Conclusions 

In this paper, we have collected all the ingredients that are needed to construct a two-loop renormalization 
procedure for the minimal Standard Model and to work out the details for constructing theoretical predictions 
for physical (pseudo-)observables. 

The procedure covers broadly (and not relying on other sources) all aspects of the problem, from diagram 
generation, to definition of the renormalization scheme, to evaluation (mostly numerically) of Green functions, 
including all customary approximations that are usually made as the small fermion-mass approximation (with 
particular attention to the analytical extraction of collinear logarithms). Our scheme can deal with arbitrary 
scales and goes beyond the introduction of MS couplings in the electroweak sector. 

Among the technical aspects which we consider to be innovative and that we have introduced in this 
paper we point out a universal ultraviolet decomposition of one- and two-loop diagrams which is closed with 
respect to multiplication and a non-negligible simplification in the structure of renormalized Green functions 
due to the (possible) introduction of a non-minimal MS prescription. We also present residues of ultraviolet 
poles, up to two-loop vertices, in terms of familiar one-loop functions. 

Techniques for expanding one- and two-loop diagrams in different regions, small fermion masses or small 
momentum transfer or both, have been introduced and discussed as well as their relevance to renormalization. 

Applications, as well as the introduction of unstable particles in the renormalization procedure, will be 
presented in a third, forthcoming, paper. 

In several appendices we have presented examples of our computational techniques. Although there is no 
substitute for the original FORM codes in modern technology, we are confident enough that the examples 
shown here suffice in illustrating our methods. The parallel development of QraphShot and of CoopBack is 
aimed to give to most complete documentation for obtaining all results which are too long to be presented 
here; if we ask the question 'might the theory in general be a lot simpler than we thought?' we may conclude 
that we have a 'nearly nifty' formulation □ showing that the Standard Model is canonically as simple as 
QED, even at the two-loop level. 
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Figure 14: The irreducible contributions to the Green function with one external photon and one external Z boson 
at O (gr 4 ) (denoted by the superscript 2). Thirteen families of Feynman diagrams give a contribution (each gets a 
subscript-label). Dotted vertices depend on /3i, Ti or the counterterms at O (<? 2 )- The black-square vertex represents 
a vertex which depends on F2 or the counterterms at O (gr 4 ) . 
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Figure 15: The correct recipe for renormalization of mass- dependent UV poles in the charged sector. 
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A Reduction of vacuum bubbles 



Following the definitions of Section 13.31 we give few examples of relations between two-loop vacuum 
bubbles and the master bubble, which we select to be Tm- For instance, we have 
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Vacuum bubbles with arbitray powers can always be reduced to a combination of products of Aq functions 
and of Ti 2 i. This results are relevant for processes where all external scales can be put to zero. 

An efficient algorithmic solution of the IBP equations for vacuum bubbles must be able to isolate special 
cases; an example is given by the following two cases: 
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We can now compute the master vacuum bubble of Eq. (|30[) . which is representable as 
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with X — x (1 — x); special cases are given in Appendix [Cl 
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B Singularities whose origin is infrared 

As mentioned in Section [2] after reduction we may have a migration of poles, from the IR to the UV; in 
this appendix we collect two examples. 
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B.l Evaluation of s°(s = ; mi,rri2,0, 0) 

In this appendix we consider the following integral, belonging to the C-family (Sect. 5 of Ref. |12j). 
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When s = and m3 = 7774 — we have V c — (1 — y) to 2 and the first term in square brackets in Eq. (|225p is 
singular at y = 1. Furthermore, 
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where we have naively assumed the relation e uv = —e 1R . The last x integration is finally performed by using 
the following relation, valid for any monomial V(x), 
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Note that the residue of the double pole has changed sign with respect to the non-IR configuration. 

B.2 Evaluation of s E (s, mi, TO2, 0, 7714, 0) 

Another example is represented by s E (s, mi, mi, 0, 7714, 0) where we have to compute the following three- 
fold integral, 
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where the quadratic form V E (see Sect. 6 of Ref. [12]) is defined as, 
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with to 2 , given in Eq. (|226p . After changing variable, y' = 1 — y, in Eq. (|230[) we use a functional relation valid 
for quadratic forms, 
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After inserting this result in Eq. (|230|) and after integration by parts we obtain the result, 
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where, once more, a = 1 — m|/s; this function, for a non-IR configuration (777,3 7^ 0), shows only a single 
pole at n = 4; in the case 7773 = (coming from reduction of integrals with irreducible scalar products in the 
numerator) there is, instead, a double pole. 
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Note that all these features arise from the reduction of S® (1, 1, • • •) and of S^(l, 1, • • •); in the case of 
pure QED all these terms are absent and, as a matter of fact, they are also absent in the doubly contracted 
Ward identity for the AA transition, even in the full standard model. It is only in the evaluation of IT(s) - 
vacuum polarization in the full standard model - that they are needed; here, for the first time, the explicit 
form of the two-loop, two-point, form factors is actually needed and represents a severe test for the final 
result. 

There are virtues and vices in any reduction algorithm; note, however, that a check on any calculation 
cannot do without reduction, i.e. without the possibility of expressing the result in terms of a basis. 



C The limit m f -> 

All diagrams contributing to a certain Green function are generated by our code QraphShot; the output 
of QraphShot is fully general and fermion masses are kept in the final answer. However, in most cases, we 
are interested in taking the limit mf — > 0, / =/= t, extracting at the same time all the collinear logarithms. 

After renormalization, we aim to compute physical observables with an algebraic-numerical approach; all 
contributions that are potentially large, e.g. collinear logarithms, are extracted with the use of an analytical 
calculation and only the collinear-free reminder is subject to a numerical evaluation. In few cases even 
this last step can be avoided given that also the finite part (in the collinear limit) is expressed in terms of 
polylogarithms. In this way we are always able to produce and explicit proof of the expected cancellations: 
no power behavior in a gauge theory etc. 

Collinear logarithms has two different origins; for instance they are present in the vacuum polarization 
function due to QED effects; for two loop irreducible contributions, however, they are sub-leading; e.g. the 
correction factor is proportional to <? 4 lnm 2 and not to <? 4 Yd 2 mf. Thus, cancellation of the leading part is 
an important test of the calculation. 

Otherwise, collinear logarithms (sometimes even power behavior) may be induced by the reduction pro- 
cedure in different pieces of the final answer, although the total is collinear free. In this cases, where there 
is no enhancement, one could as well set mf to zero from the very beginning and control cancellation of 
collinear poles in n — 4. The equivalence of the two formulations has not yet received a satisfactory answer. 

To discuss the collinear limit we introduce scaled variables, 
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and collect our results. 
- Miscellanea 



The following expansions (see Eq. (|2ip for the corresponding definitions), which have been implemented in 
QraphShot, are very useful in deriving the limit of Green functions for small values of the light fermion 
masses: 
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while, for a massive top quark, we obtain 
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- O (e) one-loop functions 
O (e) one-loop functions, Eqs. (fTl)|) - (fTT)) . are expanded as follows: 
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Vacuum bubbles 



For vacuum bubbles, see the definition in Eq. (|31[) . we introduce variables 
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and derive the following results, 
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Here, / is a generic light fermion, f ^ t. For the £ — b doublet we obtain instead, 
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Vertices and self-energies 



For all the diagrams at two loops where we have large values of one of the external \p 2 1 and small values for 
some of the internal masses it will be convenient to extract the large logarithms by defining regular functions. 
Note that individual scalar contributions may even show power behavior: such is the case, for instance, for 
Vq" when p\ 2 = and m 2 = . . . = m 5 = m when m 2 <C \P 2 \- In this case the behavior is \P 2 \/m 2 . 

In computing the vacuum polarization function we have introduced few objects that are special in the 
sense that some of the masses that appear in the argument list are zero as a consequence of the algorithm 
of reduction in sub-loops. This special feature modifies the singular structure of the function and introduces 
finite reminders that should not be confused with the regular ones for which a regulator exists; something 
that is usually mispresented is that m — 0(e regulated) is different from m — > 0. Introducing the collincar 
factor 

1 



Cl = -L f + L F L f + ^L 2 , 

we derive the following results, where {0}„ stands for n zero masses: 

s E ; fin ( S , M z ,m f , m f ,m f , m f ) = s E ; rcg (s, M z , {0} 4 ) + C f z , 
s E ' fin ( S , M zl m f , m f ,0, m f ) = s E ' rcg (s, M z , {0} 4 ) + & z , 
s E • fin ( S , 0, M w ,m f ,m f , m f ) = s E ; rcg (s, M w , {0} 4 ) + C&, 
s E ' fin ( s , M w , 0, m f ,m f , mj) = s E ; rcg (s, M w , {0} 4 ) + C&, 
s E • fin (s, M w , 0, m f , 0, m f ) = s E ; rcg ( S , M w , {0} 4 ) + C&, 



,e ; spec 



E ; speci?. 



( S ,{0} 3 ,M W ,0) + 



L f + L 2 



(243) 



- 1 



(244) 



-.e ; fi 



n (s, M w , m u , m d , m u , m d ) 



E ; rog 



(s,M w ,{0} i ) + Ci 



s E • fin (s, M w , m d , m u , m d , m u ) = s E ■ ro s( S , M w , {0} 4 ) + C£, 
s E ' fin (s, M wl m u , m d , m d , m d ) = s E ; rog (s, M w , {0} 4 ) + 

(s, M w ,m d , m u , m u , m u ) = s E '• rog (s, M w , {0} 4 ) + C™, 



e : fin 



E : fin 



(s,M z ,m u ,m u ,m d ,m u ) = s E ; rog (s, M z , {0} 4 ) + C™, 
; Kn ( S , M z , m M , m„, m„, m u ) - s E ; rcg (s, M z , {0} 4 ) + C u z , 
s E i fin (s, M z , m d , mj, m„, m d ) = s E ' rcg ( S , M z , {0} 4 ) + C d z , 
s E • fin ( s , M z , m d , m d , m d , m d ) = s E ; rcg (s, M z , {0} 4 ) + C d z , 



_,E : spec 



(s,m U) m U) 0,M W) 0) 



_,E ; specR 



(s,{0} 3 ,M w ,0) 



S * • s p° c (s, m d , m d , 0, M w , 0) = s E ; s p och ( S , {0} 3 , M w , 0) + 



3> w 



- i 



(245) 



s E;iin (s,M w ,m t ,m b ,m t ,m b ) 



; {s,M w ,m t ,0,m t ,0) 



+ 



zt-l 



L b 



x t 



— (x t L t 



,fi ; fin 



(s,M w ,m t ,m b ,m b ,m b ) 



S E;rcg ( S ,M w ,m t ,{0} 3 )-£ fa 



- lie] 

Iw • 



E : fin 



(s,M z ,m bl m b ,m b ,m b ) = s 



E;fin 



(s,M z ,m b ,m b ,m t ,m b ) 



E;rcg ( S ,M z ,{0} 4 )-L b (l-L z - l -L b 
g ( S ,M z ,{0} 2 ,m t ,0) 



«t-l 



Li. 1 - L, 



+ 



Lh Lt 



Zt 



s E • s p cc ( S , m b ,m b , 0, M w , 0) = s E > S ? CCR ( S , {0} 3 , M w , 0) + 



(246) 
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D The limit p 2 







There are several places where one needs Green functions, self-energies in particular, evaluated at p 2 = 0; 
e.g. renormalization of the electric charge requires the vacuum polarization function subtracted at p 2 = 
and corrections to G F require the W — W self-energy in the same limit. 

Of course, evaluation of two-loop two-point functions at p 2 = follows a slightly different and somehow 
independent approach. Every two-loop two-point diagram at p 2 = is a vacuum bubble that can be reduced 
to one master integral by means of integration-by-part identities. However, we have opted for a different 
approach where we start from the result at arbitrary p 2 and perform the limit p 2 — > 0; in this way we have 
a powerful check on the correctness of the original procedure. 



D.l One-loop functions 

For most of the functions appearing in an arbitrary self-energy the limit s 
nevertheless we report the corresponding expansions: 



(3 1 (z m )L (3 (z m ) 



1 

2 Z " 



1 

12 



1 

60 



with z m — m 2 /s and L m — lnz m . 

b e (2,l,s,m 2 ,m 2 ) = z' 1 
where L m = \nz m . Similarly we obtain 



Lm + - a zj (L m - 1) + — zj (L 



30 



2> 



&o(l, 1, s, m 2 , m 2 ) = -L 2 m + ^z m l L m + z m 2 (L r 



210 



is almost trivial; 

(247) 

(248) 

-) 
14 j ' 



h[ ) [l. l.*A).i l r) = -2 + 2L lll - Ll - z~ x - L w ) - z~r ( --!,„). 



,17 1 



18 



6g(2, 1, s, m 2 , 0) - 2 z m x (L m - 1) + z7 2 (L m - -). 



(249) 



The expansions for 6q with unequal masses are rather long and will not be reported here. 



D.2 Two-loop functions 

For two-loop functions we rewrite scaled functions, e.g. s A ; fin in terms of the original ones, e.g. S A ; when 
a Taylor expansion is needed we derive it in terms of vacuum bubbles. Consider an example; given 



3 A;fin (s,TOl,m2,TO 3 ) = S 1 



+ 



where the operator P uv acts as 



L - 2 ~ 3 ^ m < ) + 2Lw ^ m * Ll ~ L ™ ? m * 
Puv 5 ,A (s,mi,m 2 ,m 3 )J , 



Ptfv e _n = 0, 
we use the first few terms in the expansion, e.g. 



S ,A (s,m 1 ,m 2 ,m 3 ) = T m (mi, m 3 , m 2 ) + 7i 2 i(toi, m 3 , m 2 ) - -T m (mi,ra 3 ,m 2 ) 
+ — m3Ti3i(mi,m3,m2) s + (s 2 ) . 



(250) 
(251) 



(252) 



IBP techniques are then use to reduce vacuum bubbles to our master integral T 12 i(toi, TO3, m 2 ). Therefore, 
one function suffices in deriving the final result. This sort of results is relevant whenewer we need corrections 
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to the limit s^Oor when we want to prove that vacuum bubbles give the correct limit of the theory. For 
tensor integrals some additional work is needed; first we give an example, 



d n qi d n q 2 



{q\ + m\r [q\ + ™§)" 2 ((qi fe) 2 + m 2 )" 3 
and derive the reduction of vacuum bubble form factors in terms of scalar vacuum bubbles 



(mi,m 3 ,m 2 ) (253) 



n '^n 1 ,n 2 ,n 3 ^ ^ni-l,,ri2,n 3 m l Jni,n 2 ,n3i 

R ^iii,n 2; ti3 ~ 2 ^"i- 1 ii™2,n3 + 2 ■^"i> n 2 — l,ns — 2 ^ n b n 2i n 3- 1 + 2 H 3 _ 777-1 _ ^™ 1 ' r ' 
n ^n 1 ,n 2 ,n 3 ~ ^«i,n2- 1,713 — TO 2 ^mi ,"2 ,«3 1 

where all T have argument toi, to 2 , to 3 . Suppose now that we have to compute 



{qi + mf} ((qi+p) 2 +m\) ({q x - q 2 f + to 2 ) (<? 2 + to 2 ) ((<72+p) 2 + ™ 3 ) 



(254) 



(255) 



in the limit p 2 — > 0. the recipe is simple: expand the propagators, e.g. 



((ft +p) 2 + m 2 )^ = (g 2 + m 2 )^ - (p 2 + 2p • gi ) (qf + mj)-"- 1 + . . 



(256) 



to obtain 

nSf^O, {to}) = -T 131 (m ll m 3 ,m 2 ) - 3 T 22 i(mi, to 3 , m 2 ) 

+ (to 2 - to 2 + to 3 ) T 2 3i(toi,to 3 , to 2 ) + 2 to 2 T 32 i(toi, TO 3 ,TO 2 ) - A (2, TOl ) A (3, TO 3 ), 

n5f 2 (0, {m}) = -3T 22 i(toi,to 3 ,to 2 ) + 2 to 2 T 23 i(toi, to 3 , to 2 ) 

- r 3 ii(mi,TO3,TO 2 ) + (to 2 - to 2 + to 3 )T 32 i(toi,to 3 ,to 2 ) - A (2, to 3 ) A (3, mi). (257) 

For tensor integrals of higher rank we need to define the T^ kl form factors as 



/ 



d n qi d n q 2 



M v a 

Qi ij % ii 



(ql+mlr {4+mir ((<Zi - + m 2 2 ) 



+ T B ;ni ,n 2 ,n 3 (™4 , ™3 , ™2 ) 

+ ^^.n, Hi ,™3,m 2 )<W<W (258) 



with the following results: 



2 rpll 22 



1122 



T 1122 
"•^c; m,n 2 ,n 3 



A ; i»i ,n 2 ,ii3 

~t~-^ni— 1,712 — 1,7.3 ^2 -^Tll- 1,712)713 ^1 -^7.1 ,7.2 — 1,7.3 ~l~ ^1 ^2 -^7.1 ,7.2 ,713 ) 



(259) 



/ _ i \ rrill22 — Cl — n 2 "l 7^H22 4_ _ T 

l n !J 1 c;n 1 ,n 2 ,n 3 ~ I 1 n J J B ; rii ,n 2 ,n 3 "+" , 1 n\ -2/ 



+ (^ - ^) ^711-1,712-1,713 - ^ ^1-1,7.2,7.3-1 + (^ ~ ) ™2 + ^ (™3 ~ TO 1 ) 



[ 711— 1,712, n3 



1 ^ r ,^ l.o 9 9 \ i-n ^ 1 99 ^ \ / 9 2 2 \ 

+ -;T n n n3 _ 2 + - (to x +to 2 - to 3 )T„ „ „ 3 _i + TO 1 m 2 + - A(m 1 , m 2 , to 3 ) 

4 2 L n 4 



4r. 



711,712 — 2,713 



It 



711,712 — 1,713 — 1 



+ 



/l 1 \ 2,1/2 2\ 

(- - - mj + - m 3 - to 2 ) 
n 2 2 



T 



7.1,712 — 1,713 



(260) 
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n(n- l)(n + 2)Tj 



1122 

B ; ni,n 2 ,n 3 



• flT, 



2 "^ni- 1,712, "3 — 1 



Tii— 2,n 2 ,"3 

2 , 1 / 2 



(1 



n) T, 



Til — 1,712 — 1,T13 



- !) TO 2 + ^ ( m l - TO 3) r ni-l,n2,n 3 + ^ nT ni ,n 2 -2,n 3 - ^ nT ni;n2 _i ;Tl3 _i 



(^n-l)m? + ^(TO2-m§) 
(1 — n) m\m\ — — n \(m\, m\, m\) 



_t_ ^ rp 1/2 2 2\rrn 

-l,7l 3 ' ~£ 12 1 711,712,713 — 1 ~ ^ 71 \ m 3 ~ m l _ m 2 J J Til ,712 ,713 — 1 



T 



(261) 



2 T 1222 _ T 1222 _ T 1222 

11 J -A;ni,n 2 ,n 3 ^ B ; rii ,n 2 ,n 3 ,t - i c;n 1 ,n 2 ,n 3 

+ 2 -^"1— 1 !™2-l,«3 _ 2 777-2 ^™ 1_1, ™ 2 ' rl3 2 "^™1'™2-2,71 3 — 2 ^7li,7l 2 — l,7l 3 -l 

+ ^ ( m 3 - 2m 2 - TOi)r ni)n2 _i in3 + ^m|T nii „ 2i „ 3 _i + (ml + ml - mj) T„ lirt2 , n3 (262) 



nl222 



L C ; 711,712,713 ( n + 1) -^B ; Tli,n 2 ,Tl 3 + 2 ~~ 1 '™ 2 ~~ 1 >™ 3 



"q ^2-^711—1,712,713 + q ^Tll ,712 — 2,713 r, ^711,712 — 1,713 — 1 



+ - 



/ 2 2\ 2 

(m 3 - mj - m 2 



T, 



2n 



2 rri 
^2 Til ,712,713 



n lt 712 — 1,713 

2 f m 2 2 



-1 - m 2 ( TO 3 ~ m l - m 2) r ni,™ 2 ,n 3 : 



(263) 



n(n + 2)Tif n 2 i>n2i 



«3 



2 -^711 — 1,712 — 1,713 
+ 2 - 7n 7ii,n 2 -2,n 3 — 2 -^711,712—1,713 — 1 ~~ 



2 ^2 -^Tll— 1,712, "3 

1 2 , 2 1 2 

-m, + m 9 m, 

2 2 J 



2711,712— 1,713 



+ ^ ^2 7 1 n 1 ,„ 2 ,„ 3 -i - ^ m 2 ( m 3 - ml- ml) T nu 



(264) 



2 T-.2222 __ rpllll _ rpllll 

i A;n 1 .n 2 ,n 3 B \ 7i\ ,7i 2 ,713 ' fc J c ; ni ,n 2 ,n 3 

+T Tli rl2 _2,n 3 — 2 }7l2 2~7ii,7l 2 — 1,713 + ?7l2 T" ni ,n 2 ,n 3 7 (265) 



T-.2222 _/ 1 1 N T 2222 , J_ T 

J C; 711,712,713 V' "T" JJ J B ; 711 ,712,713 ' „ J «l ,n 2 -2,ri 3 

— 2 m 2 T" ni .n 2 — l.n 3 H ^2 -^711 .n 2 ,n 3 j (266) 

n n 



n (n -\- 2) T B . ni 7l2 n3 — 7rj iri2 _2,n 3 2 m 2 T ni n2 — \ n3 + m 2 T rii r j 2 ri3 . 



(267) 



For instance, we obtain 



S 2 E 2i(0,{m})^4(i--^-) 
n n + 2 



T 131 -2mlT 141 +m*T 151 



(268) 



Laurent expansion 



There are cases, however, where a Taylor expansion is not allowed. Consider S ,c (s,mi,ra 2l {0}2) which 
contains propagator factors q 2 and (q2 + p) 2 ■ The limit p 2 — > is IR divergent and the corresponding 
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singularity should not be regularized by e, a fact that would change the pole structure of the diagram, but 
a Laurent expansion around p 2 = is needed. In the following list, where we define 



M 2 
In w r 



In- 



s M 2 ' 

and where L an is given in Eq. (|2ip , we present the relevant limits 

'23 



In 



M 2 ' 



(269) 



S c ' hn (s, 0, to, {0} 2 ) = L w L an - {L w + L all ) 
+ 2C m + l -C 2 rn + 0{ S ), 



5 C;fin (s,m 1 ,m 2 ,{0} 2 



L 



"23 . x 2 

6 X\ — X2 



l -L 2 

2 w 



--C -(£ -C ) X2 

r . *-"mi \L-"m\ ^??i2 ) 
D X\ — X2 



l^+3C(2)+3£ mi 



1 



■1-2 



-C(2)-l-2T 1 fi 2 n 1 (0,m 1 ,m 2 )] + O (s) . 
Note that these functions are of O (in 2 s) for s — > 0. Furthermore, 



(270) 



S 5 E ; &n ( S , 0, TO, {0} 3 ) = L w L an - (i W - Lan) ( - - L wi 



s S E ' fin (s, m, to, {0} 3 ) = L tt . L an - (L lv - L an ) 

s S E ' fin (s, m 1 ,m 2 , {0} 3 ) = L w L an - (L w - L an ) 

1 



- ( T 2 A- T 2 



c 2 



r 
r 



W 

X2 



T2 , pi 



11 49 

6 m 12' 

6 m 4' 



xi — x 2 



« ( + ^an + ^mi ) _ To + "ft" 



49 11 



x 2 



* 



12 6 

— (c - C \ --[ r ~ - 11 



(271) 



D.3 Limit m/ — * of p 2 

After performing the limit p 2 — > we still would like to set light fermion masses to zero, whenever 
possible. Such is the situation for our main traget, the inclusion of a and of G F in the IPS. Functions for 
which a Laurent expansion is needed include all cases where we have a vacuum bubble T nin27l3 (0, to, M), or 
permutations of it, and n% > 2. For expanding these functions we have developed a special technique, based 
on the idea of writing differential equations for diagrams |24j . This technique works in all cases where we 
know how to compute integration constants. 

- T121 (to, M, to) with to -> 

Consider the following example: 7i2i(to, M, m) with to — * 0. We write the equation 

d 



dm- 



■ Ti2i(to, M, to) = — T 22 i(m, M, to) — Ti 22 (to, M, to). 



(272) 



We use IBP identities to express the r.h.s of Eq. (|272p in terms of Ti 2 i(to, M, to). Using Eq. (|3"Tj) we make 
the ansatze, 

/ 2 \ n / 2 \ n / 2 \ n 

T f (m,M,m) = J2a n (^- 2 ) + C rn £ 6„ (™ ) + ^ £ c„ |U , (273) 
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insert into Eq. (|272[) and derive 



1 



a = I>(0, M, 0) = -| - | C(2) + A, - Cl, (274) 

where C m = In m 2 /M 2 and £ M = In M 2 /M^. Here, the integration constant is known, see Eq. (|242p . 
Furthermore, we have 

1 



a i 

as 



2 + 2£„, 
1 



a 2 



-2C(2)-5-2£ M -2£t-4a 



92 



- 12 C(2) - — - 8 C M - 12 £ 2 f _ 24 a 



6o = 0, fei = -2, 
c = 0, ci = 0, 



6 2 = -3 + 2£ M , 6 3 = 
C2 = -1, c 3 = -4, 



32 
~3 



-— + 8A 



etc. 



Tf(m, M, 0) with m — > 



Similarly to the previous case, we obtain 

/ 2 \ ' 

T f (m,M,0) = J2 a " \ w) 

m — n V / 



M 2 

W n =0 



2 \ n 

m 2 ; ' 



with coefficients ao = 7/(0, M, 0) again given in Eq. (|242p and 

1 /l 
2 



ax = 1 + A*, 



«2 



A 



1 /l 



«3 



A 



60 = 0, 

small, unequal, masses 



bx = -l, b 2 = --, 



For small, unequal, masses we write two differential equations, 



(275) 



(276) 



(277) 



T 12 x{m 1 ,M,m 2 ) = -T 22 i(mi, M, m 2 ), 



dm 2 

and introduce the ansatz 



dm 2 



Ti2i(mi,M,TO 2 ) = -Ti2 2 (mi,M,m 2 ), 



(278) 



mi + !)!n 



Tf(m 1 ,M,m 2 ) = a Q + ai — j-^ 



a 2 ■ 



4 i 4 



M 4 



a 3 



M 4 



04 ■ 



m? + m 6 



A/ 6 



2 2 TO 1 + TO 2 

+ a 5 mj m 2 — — 



+ ln 



In 



M 2 



M 2 W 



2 



2"n 



&1 M^ + 62 M 4 



4 9 4 



+ b 5 
+ b 5 



ln4iln mi 



M 2 M 2 



ci 



/^l //^2 

M 4 



+ C2 m 1 m 2 ■ 



M 6 



The solution is a = 7/(0, M, 0) (Eq.(2I2J|) and 



«x = 


1 + A* 






1 


A 


04 = 


3 U 4 






61 - 


-1, 


62 


Cl = 


-1, 


c 2 



«2 



2 U" 1 "^ ' • 



1 ~T~ ^Af 7 

-2. 



4 2 

m\m 2 






M 6 




4 2 




m|i 


M 6 


M 6 . 



a 3 = 2 + C(2) + Ai + 2 ao , 
a 5 = 5 + 2C(2) + £ M +2£^+4a , 

63 = -^, 6 4 = -l + 2£ M , 6 5 = -4 + 2£ A 



(279) 



b 6 = 



(280) 
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- one small mass 
For one small mass, m <C mi, m 2 , we obtain 

— I + c m J2 bn (ttt) > ( 281 ) 

,, ,, w / n= i V / 

with a solution ao = 7/(0, mi, 777,2) and 

r r X2 _l r f X2 . 1 ^ . X2 r X2 _l /Yo^ 3:2 _l ^ _l 1 . o X2 

fli = A£2^ + A ^"x 2 " + xj + £ i ^ ~ £2 ^ + c(2) ^ x 2 x +2a °x 2 ' 

», - r X2 r X2 1 
& i- £ i^-^ 2 ^-^. 



x 3 y V 2 x 4 2 x 2 y 1 V x 4 x 3 y 

I 7 Xn Xo \ , . In , . £9, 3?9 9 I? 11 ( X% Xo 

^ 2 [-mk-lk +3C(2)^ + 2C(2)^ + 3-| + --4 + --,+ao 6-^+4 " 



2 X 4 X 3 y sv ; X 4 sw X 3 I 4 2I 3 4I 2 u V X 4 X 3 



X 4 X 3 7 V X 4 X 3 / X 3 2 X 2 



x 2 


-2^- 


X 4 


X 3 







/ ,2 
■ Xn 

4 

X 5 


+ 2 


^2 

X 4 


^2 


2:2 




X 5 


X 4 


) 
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x^ + 


~6~ 


X 4 



a3 = A£2|10 £| + 12 £| + 3£|) +£l (_||_ 4 i + 2 £l + iJ. 3 

^9 „ ^9 „ ^2 ^ Xn 59 49 s 2 11 

10C(2)^ + 12C(2)^ + 3C(2)^ + 10^ + y ^ + y ^ + -- 3 
a0 (20^ + 24|| + 6^ 

& 3 = A(10|| + I2f ^^(-lO^-ia^-S^-lO^-T^-i^, (282) 
where we have introduced 

9 9 
77Z - / / I 

aj» = tt|-, X = xi - X2, d = Inxi, C m = In (283) 



E Finite mass renormalization in Tl R 

As we have discussed in the previous sections there is full control upon finite parts at arbitrary scale. In 
the following Mg = M/cg, M being the renormalized W boson mass; furthermore, M H (etc) is the H boson 
(etc) renormalized mass. Consider, for example, the subtracted two-loop photon self-energy, as defined in 
sect. 6 of I, 

n*(s) = H« ;ext ( S ) - <) ;cxt (0) = 1 Ff» + F* A (s) + s 2 F^is). (284) 

Before using n R (s) for physical predictions all renormalized masses must be replaced by on-shell masses, 
using M R = M os + O (g 2 ) in one-loop terms and M R = M os in two- loop terms (see Section [531 Eqs. (|100[) . 
(j9l)|) and (|104[1 ). For instance, the effect of finite mass renormalization on H(s) — H R (s) — II R (0) is as follows; 
define UV finite factors 

M 2 = M£ ( 1 + ^ SM * ) ' ( 285 ) 
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etc. In the massless limit for fermions we get 



n(D( s )^n( 1 )( s ) + 4^An( s ), 



(286) 



An(s) = ^ E 8mf + ±x w (l + 12x w ) SM 2 b^(2,l, Sl M w ,M w ) 

- - 4Am?6g n (2,l, s,m t ,m t )-j Uw + j) £A/ 2 + | fx r + i ) 5m?, 



(287) 



with Xi = mf/s. The complete expression for ±I(s) is too long to be reported here and can be found at 



http : / / www . t o . infn . it / ~ giampier /REN /pis . ps 
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